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Near-infrared spectroscopy of Near-Earth Objects (NEOs) connects diagnostic spectral features to specific 
surface mineralogies. The combination of spectroscopy with albedos and diameters derived from thermal 
infrared observations can increase the scientific return beyond that of the individual datasets. For 
instance, some taxonomic classes can be separated into distinct compositional groupings with albedo 
and different mineralogies with similar albedos can be distinguished with spectroscopy. To that end, 
we have completed a spectroscopic observing campaign to complement the ExploreNEOs Warm Spitzer 
program that obtained albedos and diameters of nearly 600 NEOs (Trilling, D.E. et al. [2010], Astron. J. 
140, 770-784. http://dx.doi.Org/10.1088/0004-6256/140/3/770). The spectroscopy campaign included 
visible and near-infrared observations of ExploreNEOs targets from various observatories. Here we pres- 
ent the results of observations using the low-resolution prism mode (~0.7-2.5 pm) of the SpeX instru- 
ment on the NASA Infrared Telescope Facility (1RTF). We also include near-infrared observations of 
ExploreNEOs targets from the MIT-UH-IRTF Joint Campaign for Spectral Reconnaissance. Our dataset 
includes near-infrared spectra of 187 ExploreNEOs targets (125 observations of 92 objects from our sur- 
vey and 213 observations of 154 objects from the MIT survey). We identify a taxonomic class for each 
spectrum and use band parameter analysis to investigate the mineralogies for the S-, Q-, and V-complex 
objects. Our analysis suggests that for spectra that contain near-infrared data but lack the visible wave- 
length region, the Bus-DeMeo system misidentifies some S-types as Q-types. We find no correlation 
between spectral band parameters and ExploreNEOs albedos and diameters. We investigate the correla- 
tions of phase angle with Band Area Ratio and near-infrared spectral slope. We find slightly negative Band 
Area Ratio (BAR) correlations with phase angle for Eros and Ivar, but a positive BAR correlation with phase 
angle for Ganymed. The results of our phase angle study are consistent with those of (Sanchez, J.A., Reddy, 
V., Nathues, A., Cloutis, E.A., Mann, P., Hiesinger, H. [2012], Icarus 220, 36-50. http://dx.doi.org/10.1016/ 
j.icarus.2012. 04.008, arXiv:1205.0248). We find evidence for spectral phase reddening for Eros, Ganymed, 
and Ivar. We identify the likely ordinary chondrite type analog for an appropriate subset of our sample. 
Our resulting proportions of H, L, and LL ordinary chondrites differ from those calculated for meteorite 
falls and in previous studies of ordinary chondrite-like NEOs. 

© 2013 Elsevier Inc. All rights reserved. 


1. Introduction 

Since many asteroids have diagnostic spectral features in near- 
infrared wavelengths, we can investigate their mineralogies remo- 
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tely. These connections to mineralogy allow for detailed investiga- 
tions of individual objects and can be applied to questions 
regarding the compositional distributions of asteroids throughout 
the Solar System. In near-Earth space, one puzzle is the relative dis- 
tribution of ordinary chondrite-like material. Previous work (e.g. 
Vernazza et al„ 2008; Dunn et al., 2013; de Leon et al„ 2010) sug- 
gests that a majority of Near-Earth Objects (NEOs) have composi- 
tions similar to LL ordinary chondrites. This result was 
unexpected since LL ordinary chondrites are only ~7.9% of all 
meteorite falls (Burbine et al., 2002). The gap between these two 
populations will likely be bridged as we start to observe smaller 
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objects, the direct precursors to our meteorite collection, and as 
our spectroscopic observing sample sizes grow. 

There are several quantitative approaches to near-infrared 
spectral analysis. Taxonomic classification based on measurements 
of the spectral slope and absorption features is often used as a first 
step. Basic compositional can be inferred if the average composi- 
tions of each of the taxonomic classes are known. Taxonomy can 
suggest the presence of certain minerals, such as mafic silicates, 
but in many cases the mineralogical interpretation can be ambigu- 
ous. For those objects with broad 1-pm and 2-pm absorption fea- 
tures we can study the composition more thoroughly using band 
parameter analysis. Band parameter analysis has been used in 
many studies (e.g. Gaffey et al., 1993; Thomas and Binzel, 2010; 
de Leon et al., 2010; Dunn et al., 2013; Burbine et al., 2009; Mosko- 
vitz et al., 2010; Reddy et al., 2011a) to determine the compositions 
and relative abundances of the olivines and pyroxenes present on 
the surface of various asteroids. For olivine-pyroxene mixtures, 
the band center wavelength of the broad 1-pm absorption feature 
is a function of the relative abundance and composition of olivine 
and pyroxene, and the band center of the broad 2-pm absorption is 
a function of the pyroxene composition. Additionally, the Band 
Area Ratio (the ratio of the area of the 2-pm band to the area of 
the 1-pm band) is a function of the relative olivine and pyroxene 
abundances. Using the derived compositions, the diversity and dis- 
tribution of material in near-Earth space can be investigated 
remotely. 

The regions within the parameter space defined by calculated 
band centers and Band Area Ratios provide insight into meteorite 
analogs and the diversity of compositions among the asteroid sam- 
ple. These regions are correlated to potential meteorite types using 
laboratory derived mineralogies of various meteorite analogs. The 
relationship between band parameters and mineralogical composi- 
tion is particularly well studied for ordinary chondrites and basal- 
tic achondrites (e.g. Dunn et al., 2010; Burbine et al., 2009, 2007; 
Reddy et al., 201 1 a, b, 2012b). 

The combination of spectroscopy and albedos and diameters 
derived from thermal infrared observations can increase the scien- 
tific return beyond that of the individual datasets. To that end, we 
have completed a ground-based observing campaign to measure 
visible and near-infrared spectra of near-Earth objects in support 
of the Warm Spitzer project entitled “The Warm Spitzer NEO Sur- 
vey: Exploring the history of the inner Solar System and near Earth 
space” or ExploreNEOs (Trilling et al., 2010 and Trilling et al., sub- 
mitted for publication). This two year program (2009-2011) was 
allocated 500 h to determine the albedos and diameters of nearly 
600 NEOs using the Warm capabilities of the IRAC (Infrared Array 
Camera, Fazio et al., 2004; Werner et al., 2004). In this paper, we 
present our near-infrared observations from the SpeX instrument 
at the NASA Infrared Telescope Facility (IRTF). We supplement 
our observations with spectra of ExploreNEOs targets from the 
MIT-UH-IRTF Joint Campaign for NEO Spectral Reconnaissance. 


2. Data 

2.3. Observations and reduction 

The observations presented in this paper were obtained using 
the low-resolution prism mode of the SpeX instrument (Rayner 
et al., 2003) on the NASA Infrared Telescope Facility (IRTF) between 
1 September 2009 and 28 January 2012. Each object was observed 
with the dichroic out of the light path, which resulted in spectra 
covering ~0.7-2.5 pm. We used the 0.8 x 15 arcsec slit with the 
abba 7.5 arcsec dither pattern. Observations were taken at the par- 
allactic angle. Each individual exposure was nominally limited to 
120s due to variability of the atmosphere. Brighter object exposure 


times were calculated to avoid any non-linearity in the detector. 
Solar standards were selected using a list provided by S.J. Bus (per- 
sonal communication) and SIMBAD searches for G-dwarfs with so- 
lar-like B-V and V-I< colors near the position of the object. We used 
the internal SpeX calibration macro for flat-field frames and wave- 
length calibration argon lamp frames. 

The spectra were reduced in one of two methods that employed 
standard near-infrared reduction techniques. One method used 1DL 
tools (e.g. Emery et al., 2011) to perform the reduction, while the 
other used the IRAF apextract package with telluric atmosphere 
correction by the ATRAN model atmosphere (e.g. Lord, 1992; Riv- 
kin et al., 2004). Both techniques included the creation of bad pixel 
maps, division by the flat field, subtraction of the dithered pairs to 
remove atmospheric emission, wavelength calibration, and the 
extraction of two-dimensional spectra to one-dimensional spectra. 
We chose solar-like stars as telluric calibrators in order to also cor- 
rect for the spectroscopic signature of reflected sunlight from the 
targets. These techniques were tested against each other on a small 
number of objects and the reductions were nearly identical. Some 
reduced spectra show evidence of telluric absorption in the ~1.4 
and ~1.8-1.9 pm regions, which should not be confused with fea- 
tures that present compositional information. The one exception to 
these described reduction techniques is the analysis of (3552) Don 
Quixote. Mommert et al. (submitted for publication) discovered 
cometary activity on Don Quixote at the time of our ExploreNEOs 
observation. Therefore, a small reduction artifact found in the 
spectrum needed to be addressed. The final spectrum presented 
in this paper and in Mommert et al. (submitted for publication) 
was reduced using Spextool (Cushing et al., 2004) and the ATRAN 
atmospheric correction. 

We present 125 observations of 92 ExploreNEOs targets (see Ta- 
ble 1 and Fig. 9). Some of these spectra were used to taxonomically 
classify objects for the study of average albedo by taxonomic com- 
plex presented in Thomas et al. (201 1 ). Additional visible and near- 
infrared spectra from the campaign will be presented in future 
papers. 

2.2. MIT-UH-IRTF Joint Campaign for NEO Spectral Reconnaissance 

A number of ExploreNEOs targets were also observed with SpeX 
as part of the MIT-UH-IRTF Joint Campaign for NEO Spectral Recon- 
naissance. 2 The objects were identified using the MIT catalog of 
asteroid spectra and include observations from the MIT Joint Cam- 
paign, Binzel et al. (2001), Binzel et al. (2004a), and Rivkin et al. 
(2004). These objects were also observed with the low-resolution 
prism mode on SpeX with the 0.8 x 15 arcsec slit and the 7.5 arcsec 
dither pattern. These objects were reduced using the IRAF apex- 
tract package and the ATRAN model atmosphere. Data taken prior 
to mid-2007 (prior to and including sp6l) used the dichroic to re- 
move light short ward 0.82 pm. The blue end of the data for subse- 
quent observations ranges between ~ 0.7-0.8 pm. Visible 
wavelength data from the Small Main-Belt Asteroid Spectroscopic 
Survey (SMASS) are included when available (Xu et al., 1995; Bus 
and Binzel, 2002). 

We include 213 observations of 154 ExploreNEOs targets. Ta- 
ble 2 shows all the spectra used in this study, including the file 
name from the MIT-UH-IRTF Joint Campaign. The file names are 
marked with a footnote if the file does not include data at or below 
0.75 pm, which affects our band parameter analysis (Section 2.4). 

The combination of our survey and the MIT-UH-IRTF Joint Cam- 
paign for NEO Spectral Reconnaissance (hereafter “MIT Joint Cam- 
paign”) yields a total of 340 observations of 187 ExploreNEOs 
targets. 


2 http://smass.mit.edu/minus.html. 
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Table 1 

Table of observations of ExploreNEOs targets from our spectroscopic survey. The table includes object number, name, the UT date and time of observation, the integration time 
(Tint) in minutes, airmass, the solar standard used in the reduction, and the taxonomic classification. 


Number 

Name 

UT 

Tint 

Airmass 

Standard 

Taxonomic type 

433 

Eros 

01 -September-2009 12:06 

2 

1.20 

SAO 127422 

Sq/Q 

433 

Eros 

01 -November-2009 05:42 

18 

1.05 

LI 13-276 

Q/Sq 

433 

Eros 

06-Mar-2010 04:52 

8 

1.32 

SAO 75438 

Q/Sq 

433 

Eros 

28-January-2012 12:17 

1.07 

1.08 

HD 91688 

Q/Sq 

1036 

Ganymed 

19-Mar-2010 09:50 

16 

1.42 

HD 111017 

Sr/S 

1036 

Ganymed 

10-June-201 1 11:48 

2 

1.14 

HD 334883 

Sr 

1036 

Ganymed 

04-July-2011 10:31 

2.7 

1.28 

HD 204570 

Sr 

1036 

Ganymed 

17-September-2011 14:08 

3 

1.46 

Hya 64 

S 

1036 

Ganymed 

26-September-2011 10:21 

2.7 

1.24 

HD 232456 

Sr 

1036 

Ganymed 

19-October-2011 07:48 

2.7 

1.28 

SAO 75105 

Sr/S 

1627 

Ivar 

1 9-February-201 0 10:26 

32 

1.01 

HD 89547 

s 

1864 

Daedalus 

07-Mar-2010 14:41 

32 

1.42 

HD 132076 

Sq/Q 

1865 

Cerberus 

05-August-2010 09:08 

96 

1.08 

Ll 10-361 and LI 12-1333 

s 

1866 

Sisyphus 

01 -May-2011 11:56 

16 

1.05 

HD 144684 

s 

1866 

Sisyphus A 

10-June-2011 06:16 

16 

1.18 

HD 125982 

s 

1866 

Sisyphus B 

10-June-201 1 06:52 

16 

1.13 

HD 125982 

s 

1866 

Sisyphus C 

10-June-201 1 07:19 

16 

1.11 

HD 125982 

s 

1866 

Sisyphus D 

10-June-201 1 08:22 

16 

1.14 

HD 125982 

s 

1866 

Sisyphus E 

10-June-2011 09:07 

8 

1.21 

HD 125982 

s 

1866 

Sisyphus F 

10-June-201 1 09:27 

16 

1.30 

HD 125982 

s 

1866 

Sisyphus A 

24-June-201 1 06:08 

16 

1.16 

HD 123765 

s 

1866 

Sisyphus B 

24-June-201 1 06:40 

16 

1.17 

HD 123765 

s 

1866 

Sisyphus C 

24-June-201 1 07:04 

16 

1.19 

HD 123765 

s 

1917 

Cuyo A 

24-June-201 1 10:16 

16 

1.01 

HD 342904 

Sv 

1917 

Cuyo B 

24-June-201 1 10:46 

24 

1.01 

HD 342904 

Sr 

1917 

Cuyo C 

24-June-201 1 11:31 

24 

1.02 

HD 342904 

Sv/Sr 

1917 

Cuyo D 

24-June-201 1 12:11 

24 

1.05 

HD 342904 

Sr 

1917 

Cuyo E 

24-June-201 1 12:48 

24 

1.11 

HD 342904 

Sv/Sr 

1917 

Cuyo F 

24-June-201 1 13:29 

24 

1.22 

HD 342904 

Sr 

1917 

Cuyo 

14-August-2011 06:59 

48 

1.01 

Ll 07-684 

Sr 

1943 

Anteros 

02-September-2009 11:30 

32 

1.15 

SAO 107874 

Sq 

2102 

Tantalus 

04-July-2011 06:36 

16 

1.17 

HD 234075 

Sr 

3103 

Eger 

24-June-201 1 14:38 

40 

1.06 

HD 206658 

C/X INDET 

3122 

Florence 

05-August-2010 15:26 

10 

1.03 

Ll 15-271 

Sq/Q 

3200 

Phaethon 

01 -September-2009 13:53 

64 

1.20 

SAO 3963 

B 

3552 

Don Quixote 

18-0ctober-2009 05:37 

20 

1.54 

Ll 13-276 

D 

3554 

Amun 

01 -February-2010 14:44 

32 

1.12 

SAO 97431 

D 

4015 

Wilson-Harrington 

02-September-2009 05:02 

48 

1.40 

HD 148116 

C* 

4055 

Magellan 

01 -September-2010 05:13 

32 

1.06 

HD 156486 

V 

4283 

Cuno 

19-October-2011 12:17 

32 

1.07 

HD 16397 

Sq/Q 

5143 

Fleracles 

27-September-201 1 12:58 

32 

1.04 

HD 275948 

Q 

5587 

1990 SB 

28-January-2012 14:25 

48 

1.25 

SAO 157011 

S/Sr 

5626 

1991 FE 

02-September-2009 06:47 

32 

1.25 

HD 190524 

Sq 

5646 

1990 TR 

01-November-2009 14:22 

48 

1.03 

HD 60913 

Q 

5786 

Talos 

05-August-2010 11:18 

72 

1.02 

Ll 10-361 and Ll 12-1333 

Sq 

6239 

Minos 

01 -September-2010 11:49 

16 

1.12 

HD 220685 

Q/Sq 

7358 

Oze 

13-August-2011 05:44 

48 

1.32 

SAO 159708 

s 

7358 

Oze 

14-August-2011 05:47 

48 

1.32 

Ll 07-684 

s 

7358 

Oze 

17-September-2011 05:02 

48 

1.39 

Ll 10-361 

s 

7822 

1991 CS 

06-Mar-2010 05:37 

48 

1.09 

HD 288180 

s* 

8566 

1996 EN 

02-September-2009 14:58 

32 

1.13 

HD 30947 

V 

8567 

1996 HW1 

26-September-201 1 15:11 

16 

1.02 

HD 247554 

Q 

11066 

Sigurd 

1 9-February-201 0 12:32 

48 

1.10 

HD 129290 

Sr 

12711 

Tukmit 

01 -September-2009 05:07 

48 

1.25 

BD + 12 2836 

Q 

12923 

Zephyr 

04-April-2010 12:19 

64 

1.19 

HD 125982 

s 

16834 

1997 WU22 

01 -September-2009 12:24 

48 

1.14 

SAO 127422 

Sq 

17274 

2000 LC16 

01 -September-2009 09:07 

16 

1.05 

HD 205027 

D 

18882 

1999 YN4 

27-December-2010 07:52 

32 

1.33 

SAO 13263 

Sr 

19764 

2000 NF5 

27-September-2010 09:32 

32 

1.04 

SAO 109095 

Sq 

20086 

1994 LW 

24-June-201 1 08:29 

44 

1.18 

HD 133600 

C/X INDET 

21088 

1992 BL2 

19-February-2010 08:06 

64 

1.20 

HD 70996 

Q 

23187 

2000 PN9 

01 -September-2010 09:04 

80 

1.04 

HD 206658 

Sq 

24761 

Ahau 

06-January-2010 15:07 

16 

1.20 

HD 96995 

C/X INDET 

27346 

2000 DN8 

01 -February-2010 08:50 

48 

1.25 

HD 43713 

Q 

36284 

2000 DM8 

02-February-201 1 12:26 

32 

1.03 

HD 249566 and HD 89179 

Sq 

39572 

1993 DQ1 

17-September-2011 09:40 

80 

1.30 

Ll 13-276 

Q 

52750 

1998 KK17 

19-Mar-2010 14:54 

32 

1.27 

HD 175470 

V 

53789 

2000 EDI 04 

27-September-2010 05:22 

80 

1.14 

SA 115-271 and SAO 86827 

Sq 

54789 

2001 MZ7 

01 -February-2010 10:51 

16 

1.05 

SAO 97431 

C/X INDET 

54789 

2001 MZ7 

1 9-February-201 0 07:21 

16 

1.11 

SAO 60533 

C/X INDET 

65679 

1989 UQ 

15-0ctober-2010 08:19 

32 

1.13 

HD 10488 

C/Cb 

68216 

2001 CV26 

19-February-2010 14:03 

56 

1.01 

HD 129290 

Sq 

68216 

2001 CV26 

07-Mar-2010 11:47 

32 

1.06 

HD 123760 

Sq 


(continued on next page ) 
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Table 1 ( continued ) 


Number 

Name 

UT 

Tint 

Airmass 

Standard 

Taxonomic type 

68216 

2001 CV26 

19-Mar-2010 10:55 

32 

1.03 

HD 115692 

S/Sq 

68350 

2001 MK3 

02 -February-2011 05:25 

72 

1.15 

HD 23925 

S/Sr 

85839 

1998 Y04 

1 9-February-201 0 11:32 

32 

1.02 

HD 96657 

Q 

86067 

1999 RM28 

01 -February-2010 13:22 

32 

1.11 

HD 87776 

Q 

88254 

2001 FM129 

04-April-2010 09:22 

32 

1.36 

HD 82410 

a 

96590 

1998 XB 

15-December-2010 10:16 

24 

1.09 

Hya 64 

Sq 

100926 

1998 MQ 

15-0ctober-2010 11:25 

24 

1.40 

SAO 40380 

Q 

136617 

1994 CC 

18-0ctober-2009 11:49 

80 

1.12 

L93-101 and L98-978 

Sa 

137032 

1998 UOl 

15-0ctober-2010 05:06 

32 

1.02 

SAO 89822 

Sq 

137062 

1998 WM 

15-0ctober-2010 09:23 

64 

1.20 

HD 12264 

Sq 

137084 

1998 XS16 

02 -January-201 1 08:09 

48 

1.05 

HD 42160 

s 

137084 

1998 XS16 

03-January-2011 07:10 

32 

1.16 

HD 42160 

s 

137125 

1999 CT3 

02 -February-2011 13:24 

84 

1.25 

HD 92670 

Q 

137170 

1999 HF1 

01 -May-2011 12:32 

16 

1.23 

HD 177780 

C/X 1NDET 

137671 

1999 XP35 

28-December-201 1 10:57 

48 

1.16 

SA 98-978 

s 

138883 

2000 YL29 

02-September-2009 14:03 

32 

1.06 

SAO 55805 

a 

138911 

2001 AE2 

28-January-2012 09:26 

92 

1.05 

SAO 97431 

S 

141498 

2002 EZ16 

05-December-2010 09:40 

32 

1.24 

HD 30518 

Sq 

141498 

2002 EZ16 

12-December-2010 06:55 

68 

1.08 

L93-101 

Sq 

143381 

2003 BC21 

27-September-2010 08:00 

48 

1.05 

SAO 127304 

Sr/S 

152558 

1990 SA 

01 -September-2010 11:19 

16 

1.02 

SAO 90986 

Sq/Q 

152563 

1992 BF 

02-February-201 1 07:29 

72 

1.09 

HD 249566 

K/C/X/ 1NDET 

152931 

2000 EA107 

07-Mar-2010 13:13 

48 

1.22 

HD 144423 

q 

153814 

2001 WN5 

15-0ctober-2010 06:24 

48 

1.11 

HD 200565 

K/L/Sq 1NDET 

153842 

2001 XT30 

05-December-2010 10:46 

80 

1.07 

HD 288029 

C/X 1NDET 

154029 

2002 CY46 

01 -September-2010 12:23 

92 

1.01 

HD 14786 

s 

159402 

1999 AP10 

01 -September-2009 09:49 

16 

1.32 

HD 216298 

Sq 

162998 

2001 SKI 62 

1 9-February-201 0 05:27 

52 

1.10 

HD 285295 

D 

163132 

2002 CU1 1 

07-Mar-2010 10:23 

48 

1.29 

HD 77066 

C/X 1NDET 

163243 

2002 FB3 

04-April-2010 05:10 

64 

1.18 

HD 296688 

Q 

163697 

2003 EF54 

02-September-2009 10:01 

48 

1.13 

SAO 107874 

Q 

164121 

2003 YT1 

18-0ctober-2009 13:55 

28 

1.09 

L93-101 and L98-978 

V 

198856 

2005 LR3 

12-December-2010 04:50 

80 

1.15 

LI 15-271 

Sq 

198856 

2005 LR3 

16-December-2010 05:58 

48 

1.11 

LI 15-271 

q 

207945 

1991 JW 

15-December-2010 05:26 

72 

1.17 

LI 15-271 

q 

207945 

1991 JW 

16-December-2010 04:34 

64 

1.06 

LI 15-271 

a 

214088 

2004 JN13 

06-January-2010 04:51 

48 

1.34 

LI 15-271 

Sq 

218863 

2006 W0127 

18-0ctober-2009 10:23 

64 

1.14 

L93-101 and L98-978 

Sq 

219071 

1997 US9 

05-December-2010 07:41 

64 

1.34 

HD 16472 

Q 

220124 

2002 TE66 A 

04-April-2010 10:44 

8 

1.03 

HD 124019 

Sq/Q 

220124 

2002 TE66 B 

04-April-2010 14:17 

16 

1.30 

HD 124019 

Sq 

247517 

2002 QY6 

05-August-2010 06:25 

112 

1.18 

Ll 10-361 and L112-1333 

Sr 

265187 

2003 YS117 

02 -February-2011 09:58 

88 

1.15 

HD 68972 

Q 

275792 

2001 QH142 

27-January-2012 10:49 

136 

1.01 

L102-1081 

C/X 1NDET 

297418 

2000 SP43 

17-September-2011 14:47 

36 

1.04 

Hya 64 

V 

347813 

2002 NP1 

02-September-2009 08:09 

48 

1.03 

HD 199268 

Q 


2000 CO101 

04-April-2010 07:20 

64 

1.19 

SAO 117128 

C/X/T 1NDET 


2000 TJ1 

27-September-2010 10:38 

82 

1.03 

SA 115-271 and HD 5294 

Sq 


2001 HC 

04-April-2010 14:55 

32 

1.07 

HD 125982 

Sq 


2002 OS4 

02-January-2011 06:07 

64 

1.21 

SAO 39637 

Sr 


2002 RQ25 

02-September-2009 12:35 

48 

1.36 

HD 5078 

C/X 1NDET 


2005 MC 

06-January-2010 14:18 

16 

1.25 

HD 73668 

Sr 


2.3. Taxonomic classification 

Taxonomic classifications of all 340 observations from both our 
ExploreNEOs survey and the MIT Joint Campaign determined from 
the Bus-DeMeo online taxonomic classifier 3 (DeMeo et al., 2009) 
are presented in Tables 1 and 2. We use the definitions of the com- 
plexes adopted by Thomas et al. (2011) and defined in Binzel et al. 
(2004c). 

The online classification system only returns a definitive classi- 
fication when a spectrum contains visible and near-infrared data. 
This is the case for a subset of the MIT Joint Campaign data. For 
all near-infrared only observations (those without visible spectra), 
the classification system returns a variety of possible subclasses. 
The additional absorption features and slope information con- 
tained in the visible wavelength region are vital to non-ambiguous 
classification, especially for the relatively featureless C- and X- 
complexes. DeMeo et al. (2009) show that the lack of visible wave- 


3 http://smass.mit.edu/cgi-bin/busdemeodass-cgi. 


length data allow certain classes to occupy the same principal 
component space. This prevents the data from being formally clas- 
sified by principal component analysis, and all possible types are 
returned by the system. Some of these objects can be classified if 
the presence or absence of weak absorption features or small spec- 
tral variations typical to the potential subclasses are taken into ac- 
count. These features are too small to be distinguished in principal 
component space and too small to affect the residuals between the 
input spectrum and the class average spectrum, but are indicative 
of certain classifications. 

For all near-infrared only spectra, we compare the object’s spec- 
trum with the average spectra of the possible subclasses. Potential 
classifications are ruled out when the shapes, slopes, and/or resid- 
uals of the object in question are not consistent with the average 
subclass spectra. For most objects, this examination reveals the 
most likely or two most likely subclass classifications for the ob- 
ject. Two subclass classifications are given only if the two subclass- 
es are within the same taxonomic complex. If several potential 
subclass classifications within a single complex are possible 
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Table 2 

Table of ExploreNEOs spectra from the M1T-UH-IRTF Joint Campaign for NEO Spectral Reconnaissance. The table includes object number, name, the file name from the MIT survey 
catalog (smass.mit.edu), the UT date of observation, and the taxonomic classification. 


Number 

Name 

File name 

UT date 

Taxonomic type 

433 

Eros 

a000433.spl5 

17-June-02 

Sw 

433 

Eros 

a000433.spl6 

17-August-02 

Sw 

433 

Eros 

a000433.spl7 

16-September-02 

Sw 

433 

Eros 

a000433.spl01 

22-August-l 1 

Sw 

433 

Eros 

a000433.spl02 

25-September-ll 

Sw 

433 

Eros 

a000433.spl03 

24-October-ll 

Sw 

433 

Eros 

a000433.spl05 

28-December-l 1 

Sqw 

1036 

Ganymed 

a001036.sp05 

28-Mar-01 

Sr 

1036 

Ganymed 

a001036.sp51 

03-June-06 

Sr 

1036 

Ganymed 

a001036.sp80 

27-April-09 

Sr 

1036 

Ganymed 

a001036.spl03 

24-October-ll 

S 

1036 

Ganymed 

a001036.spl04 

31 -October- 11 

S 

1627 

Ivar 

a001627.sp09 

13-January-02 

S 

1627 

Ivar 

a001627.sp73 

02-September-08 

Sqw 

1627 

Ivar 

a001627.sp74 

02-0ctober-08 

Sw 

1627 

Ivar 

a001627.sp76 

03-December-08 

S 

1627 

Ivar 

a001627.sp77 

27-January-09 

Sr 

1627 

Ivar 

a001627.sp88 

22-February-10 

Sw 

1685 

Toro 

a001685.spl3 

08-May-02 

S 

1685 

Toro 

a001685.spl7 a 

16-September-02 

S 

1864 

Daedalus 

a001864.sp05 

29-Mar-01 

Sq 

1865 

Cerberus 

a001865.sp75 

31 -October-08 

s 

1866 

Sisyphus 

a001866.sp56 

21 -November-06 

Sw 

1916 

Boreas 

a001916.sp07 

14-August-01 

Sw 

1917 

Cuyo 

a001917.sp72 

07-July-08 

Sv 

1943 

Anteros 

a001943.sp09 

12-January-02 

Sw 

1943 

Anteros 

a001943.spl03 

24-October-ll 

S 

1980 

Tezcatlipoca 

a001980.sp55 

25-October-06 

Sw 

2212 

Hephaistos 

a002212.sp57 a 

22-December-06 

Sq 

3102 

Krok 

a003102.sp02 

09-0ctober-00 

Sqw 

3103 

Eger 

a003103.sp06 

21 -June-01 

Xe 

3199 

Nefertiti 

a003199.sp03 

30-January-01 

NO TAX MATCH 

3199 

Nefertiti 

a003199.sp36 

03-Mar-05 

K 

3200 

Phaethon 

a003200.spl9 

27-October-02 

NO TAX MATCH 

3200 

Phaethon 

a003200.sp34 

10-December-04 

B 

3552 

Don Quixote 

a003552.sp85 

25-October-09 

D 

3554 

Amun 

a003554.sp78 a 

02-Mar-09 

C/X 1NDET 

3554 

Amun 

a003554.sp88 a 

23-February-10 

D 

3554 

Amun 

a003554.sp89 a 

16-Mar-10 

D 

3554 

Amun 

a003554.sp98 a 

05-April-ll 

D 

3671 

Dionysus 

a003671.sp91 

10-May-10 

C/X 1NDET 

3691 

Bede 

a003691.sp68 

10-Mar-08 

Cgh 

4015 

Wilson-Harrington 

a004015.sp81 a 

22-June-09 

B 

4055 

Magellan 

a004055.sp38 

11 -April-05 

V 

4183 

Cuno 

a004183.spl03 

24-October-ll 

0. 

4544 

Xanthus 

a004544.sp58 a 

21 -January-07 

Sq/Q 

4660 

Nereus 

a004660.6 

15-Mar-02 

X* 

5011 

Ptah 

a005011.sp56 a 

21 -November-06 

Q 

5131 

1990 BG 

a005131.sp67 

15-January-08 

Sq 

5131 

1990 BG 

a005131.sp87 

17-January-10 

Sqw 

5143 

Heracles 

a005143.sp55 

25-October-06 

a 

5332 

1990 DA 

a005332.sp47 a 

13-November-05 

Sr/S 

5496 

1973 NA 

a005496.sp96 a 

06-January-l 1 

C/X 1NDET 

5587 

1990 SB 

a005587.sp05 

28-Mar-01 

Sr 

5604 

1992 FE 

a005604.sp05 

29-Mar-01 

V 

5604 

1992 FE 

a005604.sp80 

27-April-09 

V 

5620 

Jasonwheeler 

a005620.sp81 a 

22-June-09 

I< 

5626 

1991 FE 

a005626.sp84 

20-September-09 

S 

5645 

1990 SP 

a005645.spl2 a 

14-April-02 

C/X/T 1NDET 

5653 

Camarillo 

a005653.sp35 a 

08-January-05 

S/Sr 

5693 

1993 EA 

a005693.sp58 a 

21 -January-07 

Sq 

5786 

Talos 

a005786.sp08 a 

23-August-01 

Sq/Q 

6047 

1991 TB1 

a006047.sp45 

08-0ctober-05 

s 

6239 

Minos 

a006239.sp93 

07-September-10 

Sq 

6455 

1992 HE 

a006455.spll 

15-Mar-02 

Srw 

6455 

1992 HE 

a006455.spl9 

28-October-02 

Srw 

6455 

1992 HE 

a006455.spl07 

21 -April-1 2 

Sqw 

7350 

1993 VA 

a007350.sp30 a 

15-September-04 

C/X 1NDET 

7350 

1993 VA 

a007350.sp39 a 

17-April-05 

C/X 1NDET 

7358 

Oze 

a007358.spl00 

06-June-l 1 

Sq 

8567 

1996 HW1 

a008567.sp44 a 

04-September-05 

Q 

8567 

1996 HW1 

a008567.sp48 a 

22-November-05 

Q 

8567 

1996 HW1 

a008567.sp73 

02-September-08 

Q 


( continued on next page ) 
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Table 2 ( continued ) 


Number 

Name 

File name 

UT date 

Taxonomic type 

8567 

1996 HW1 

a008567.sp74 J 

02-0ctober-08 

Sq 

8567 

1996 HW1 

a008567.spl04 a 

31 -October- 11 

a 

10115 

1992 SK 

aOlOl 15.sp55 a 

24-October-06 

Sq/S 

10145 

1994 CK1 

aul994ckl.spl05 a 

30-December-ll 

s 

11066 

Sigurd 

a011066.sp88 

22-February-10 

s 

11398 

1998 YP11 

a011398.sp68 

10-Mar-08 

Sr 

11885 

1990 SS 

a011885.sp98 a 

05-April-ll 

CLASSIFIER INDET 

12923 

Zephyr 

a012923.sp92 a 

1 1 -July-1 0 

S 

14402 

1991 DB 

a014402.sp78 

02-Mar-09 

Ch 

15745 

1991 PM5 

a015745.sp82 

07-August-09 

s 

16834 

1997 WU22 

a016834.sp77 

27-January-09 

s 

16834 

1997 WU22 

a016834.sp84 

20-September-09 

s 

17274 

2000 LC16 

a017274.sp02 

10-0ctober-00 

D 

17274 

2000 LC16 

a017274.sp83 

24-August-09 

CLASSIFIER INDET 1 

19764 

2000 NF5 

a019764.sp93 a 

06-September-10 

Sq 

20086 

1994 LW 

a020086.spl00 a 

07-June-l 1 

C/Cb 

20790 

2000 SE45 

a020790.sp03 

29-January-01 

s 

21088 

1992 BL2 

aul992BL2.sp03 a 

30-January-01 

a 

22753 

1998 WT 

a022753.spl9 a 

27-October-02 

S/Sq 

22753 

1998 WT 

a022753.sp36 a 

03-Mar-05 

Sq 

23183 

2000 OY21 

a023183.sp49 a 

30-January-06 

Sq/Q 

23187 

2000 PN9 

a023187.sp05 a 

28-Mar-01 

S/Sq 

23187 

2000 PN9 

a023187.sp98 a 

05-April-ll 

Sq/Q 

24761 

Ahau 

a024761.sp87 a 

17-January-10 

C/X INDET 

25143 

Itokawa 

a025143.5 

19-February-01 and 28-Mar-01 

Sqw 

25916 

2001 CP44 

a025916.sp55 a 

24-October-06 

Sq/Q 

25916 

2001 CP44 

a025916.sp93 a 

07-September-10 

Sq/Q 

31669 

1999 JT6 

a031669.spl2 a 

1 5-April-02 

Q 

35107 

1991 VH 

a035107.sp20 

27-December-02 

Sq 

35107 

1991 VH 

a035107.sp71 

11 -June-08 

Sq 

36284 

2000 DM8 

a036284.spll 

16-Mar-02 

Sq 

36284 

2000 DM8 

a036284.sp97 

07-February-l 1 

NO TAX MATCH 

39572 

1993 DQ1 

a039572.sp73 a 

02-September-08 

Q 

52387 

1993 OM7 

aul 9930M7.sp20 a 

27-December-02 

Sr 

52750 

1998 KK17 

a052750.sp30 a 

15-September-04 

V 

52768 

1998 OR2 

a052768.sp79 a 

30-Mar-09 

L 

53430 

1999 TY16 

a053430.sp76 a 

03-December-08 

Sr 

53435 

1999 VM40 

a053435.sp55 

25-October-06 

Srw 

53789 

2000 EDI 04 

au2000EDl 04.spl 7 a 

15-September-02 

S 

54690 

2001 EB 

au2001EB.sp05 

28-Mar-01 

S 

54789 

2001 MZ7 

a054789.sp21 a 

16-Mar-03 

C/X INDET 

54789 

2001 MZ7 

a054789.sp87 a 

17-January-10 

C/X INDET 

54789 

2001 MZ7 

a054789.sp89 a 

16-Mar-10 

C/X INDET 

65679 

1989 UQ 

a065679.sp94 

13-October-lO 

c 

66251 

1999 GJ2 

a066251.sp44 a 

05-September-05 

Q 

66251 

1999 GJ2 

a066251.sp45 a 

08-0ctober-05 

Q 

66251 

1999 GJ2 

a066251.sp46 a 

31 -October-05 

Q 

68216 

2001 CV26 

a068216.sp89 a 

16-Mar-10 

Sq 

68359 

2001 OZ13 

a068359.sp67 a 

16-January-08 

Sq 

69230 

Hermes 

a069230.9 a 

16-October-03 

Sq 

85709 

1998 SG36 

a085709.sp50 c1 

30-April-06 

S/Sr 

85818 

1998 XM4 

a085818.sp91 

10-May-10 

Srw 

85839 

1998 Y04 

a085839.sp89 a 

16-Mar-10 

Sq/Q 

85989 

1999 JD6 

a085989.sp92 

12-july-10 

K 

86324 

1999 WA2 

a086324.sp67 a 

16-January-08 

Sr/S 

86667 

2000 FO10 

a086667.spl07 a 

2 1 -April-1 2 

Sq 

87684 

2000 SY2 

au2000SY2.sp30 a 

15-September-04 

Q 

88254 

2001 FM129 

a088254.sp89 a 

17-Mar-10 

Sq 

89355 

2001 VS78 

au2001VS78.spll a 

16-Mar-02 

Sr 

89355 

2001 VS78 

au2001VS78.spl2 

1 5-April-02 

Sr 

90367 

2003 LC5 

au2003LC5.sp23 a 

05-July-03 

C/X/T INDET 

96590 

1998 XB 

a096590.sp96 a 

06-January-l 1 

Q/Sq 

99799 

2002 LJ3 

a099799.spl05 a 

30-December-ll 

Sq 

99907 

1989 VA 

aul989VA.spl9 

27-October-02 

Sr 

99907 

1989 VA 

a099907.sp65 a 

11 -November-07 

S/Sr 

99935 

2002 AV4 

a099935.sp67 a 

16-January-08 

B 

100926 

1998 MQ. 

al00926.sp94 

14-October-lO 

Sqw 

108519 

2001 LF 

au2001LF.sp23 a 

05-July-03 

C 

136564 

1977 VA 

aul977VA.sp44 a 

04-September-05 

C/X INDET 

137032 

1998 UOl 

aul998U01.spl9 a 

27-October-02 

S/Sr 

137032 

1998 UOl 

al37032.sp74 a 

02-0ctober-08 

Sq/Q 

137062 

1998 WM 

aul998WM.spl9 

27-October-02 

Sr 

137170 

1999 HF1 

aul999HFl.spl3 a 

07-May-02 

C/X INDET 

137170 

1999 HF1 

aul999HFl.sp38 a 

11 -April-05 

C/X INDET 

137170 

1999 HF1 

a!37170.sp98 a 

05-April-ll 

C/X INDET 
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Table 2 (continued) 


Number 

Name 

File name 

UT date 

Taxonomic type 

137671 

1999 XP35 

aul999xp35.spl05 a 

30-December-ll 

S/Sr 

137924 

2000 BD19 

au2000BD19.sp03 a 

29-January-01 

V 

138258 

2000 GD2 

au2000GD2.spll 

15-Mar-02 

Sq 

138883 

2000 YL29 

al38883.sp84 a 

20-September-09 

a 

141052 

2001 XR1 

al41052.sp77 

27-January-09 

Sq 

141432 

2002 CQ1 1 

al41432.sp99 a 

30-April-ll 

S 

143651 

2003 QO104 

al43651.sp79 a 

30-Mar-09 

Sq/Q 

143651 

2003 QO104 

al43651.sp80 

27-April-09 

a 

144900 

2004 VG64 

au2004VG64.sp46 a 

31 -October-05 

Sq 

144922 

2005 CK38 

au2005CK38.sp37 a 

08-Mar-05 

Sq 

152558 

1990 SA 

al52558.sp93 a 

07-September-10 

Sq 

152931 

2000 EA107 

al52931.sp89 

16-Mar-10 

a 

153219 

2000 YM29 

au2000YM29.sp58 a 

21 -January-07 

C/X 1NDET 

153814 

2001 WN5 

al53814.sp94 a 

12-October-lO 

NO TAX MATCH b 

154029 

2002 CY46 

al54029.sp93 a 

06-September-10 

s 

154244 

2002 KL6 

al54244.sp81 a 

21 -June-09 

Q 

154715 

2004 LB6 

al54715.sp96 a 

06-January-l 1 

Q 

159402 

1999 AP10 

al59402.sp84 a 

20-September-09 

a 

161989 

Cacus 

aul978CA.sp21 a 

16-Mar-03 

a 

161998 

1988 PA 

al61998.sp92 a 

ll-July-10 

Q 

162483 

2000 PJ5 

al62483.sp82 a 

08-August-09 

a 

162900 

2001 HG31 

al62900.sp76 a 

03-December-08 

Q 

163001 

2001 SE170 

au2001SE170.sp33 a 

24-October-04 

Sr 

163243 

2002 FB3 

al63243.dm07 a 

20-Mar-12 

Sq/Q 

163697 

2003 EF54 

al63697.sp83 a 

24-August-09 

0. 

164202 

2004 EW 

au2004EW.sp36 a 

03-Mar-05 

C/X 1NDET 

164400 

2005 GN59 

al64400.sp73 a 

02-September-08 

Q 

175706 

1996 FG3 

al75706.sp79 

30-Mar-09 

C/Ch 

175706 

1996 FG3 

al75706.sp80 

27-April-09 

B 

175706 

1996 FG3 

al75706-obsA.spl05 

30-December-ll 

C/Ch b 

175706 

1996 FG3 

al 75706-obsB.spl 05 

30-December-ll 

C/Ch b 

175706 

1996 FG3 

al75706.dm04 

1 -December-11 

NO TAX MATCH b 

184266 

2004 VW14 

au2004VW14.sp69 a 

13-April-08 

Q 

200840 

2001 XN254 

au2001XN254.spl2 

14-April-02 

s 

207945 

1991 JW 

a207945.sp80 

27-April-09 

Sq/Q 

217796 

2000 T064 

a217796.sp85 a 

25-October-09 

Sr 

217807 

2000 XK44 

a217807.sp85 a 

25-October-09 

Q 

218863 

2006 W0127 

a218863.sp89 a 

17-Mar-10 

Sq 

219071 

1997 US9 

a219071.spl04 

31 -October- 11 

Q 

235756 

2004 VC 

au2004VC.sp79 a 

30-Mar-09 

Sr/S 

237442 

1999 TA10 

a237442.sp91 a 

10-May-10 

S/Sr 

241662 

2000 K044 

a241662.sp92 a 

12-July-10 

Sq/Q 

248818 

2006 SZ217 

au2006SZ217.sp76 a 

03-December-08 

Sr/S 

256412 

2007 BT2 

au2007BT2.sp79 a 

30-Mar-09 

C/Cb 

257744 

2000 AD205 

au2000AD205.sp71 a 

11 -June-08 

I< 

260141 

2004 QT24 

au2004QT24.sp39 a 

18-April-05 

Q 

260141 

2004 QT24 

a260141.sp99 a 

30-April-ll 

Q 

267729 

2003 FC5 

a267729.sp98 a 

05-April-ll 

s c 

274138 

2008 FU6 

au2008FU6.sp98 a 

05-April-ll 

Sq 

274138 

2008 FU6 

a274138.sp99 a 

30-April-ll 

Q 

297418 

2000 SP43 

au2000SP43.sp02 a 

09-0ctober-00 

V 

297418 

2000 SP43 

au2000sp43.spl 02 a 

25-September-ll 

V 

297418 

2000 SP43 

a297418.spl03 a 

23-October-ll 

V 

301844 

1990 UA 

aul990ua.spl03 a 

24-October-ll 

Sq 

302831 

2003 FH 

au2003fh.spl04 a 

31 -October- 11 

Q 

310442 

2000 CH59 

au2000CH59.sp67 a 

16-January-08 

Sq 


2000 CO101 

au2000col01.sp84 

20-September-09 

X 


2002 G05 

au2002G05.spl2 a 

14-April-02 

Sq 


2002 LV 

au20021v.sp81 a 

22-June-09 

S/Sr 


2002 LV 

au20021v.sp82 a 

07-August-09 

Sr 


2005 JE46 

au2005JE46nl.sp47 a 

12-November-05 

NO TAX MATCH b 


2005 JE46 

au2005JE46n2.sp47 a 

13-November-05 

NO TAX MATCH b 


2005 JE46 

au2005JE46.sp48 a 

22-November-05 

C/X/T 1NDET 


2005 UL5 

au2005UL5.sp47 a 

13-November-05 

Sq 


2007 RF5 

au2007RF5.sp64 a 

02-0ctober-07 

Sq 


a File does not include data at or below 0.75 |tm. Observations were likely taken with the dichroic in the light path. 
b Objects can not be formally classified due to the presence of a thermal tail. 

c Object spectrum is very low resolution, only the complex can be identified. Possible classifications include members of the S-complex that are not S-subtypes. 


following examination, the object is classified as belonging to that 
complex and the classification is marked as the complex letter fol- 
lowed by a star (e.g„ S*). 

If the potential subclass classifications span more than one tax- 
onomic complex, then the object is labeled as being of indetermi- 


nate complex and is marked as such. This is particularly prevalent 
for objects of C- or X-complex. These complexes are often difficult 
to separate using only near-infrared data. The subtle visible wave- 
length spectral features are often crucial to making definitive 
assignments. 
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We include two labels that indicate an object has not been clas- 
sified. If the classification system finds no matches in principal 
component space it returns “Classifier Indeterminate.” If the clas- 
sifier returns several options in principal component space that 
are not good matches to the spectrum, we do not classify the object 
and label it as “No Tax Match”. These two labels appear in Tables 1 
and 2. The problems classifying the “No Tax Match” spectra are 
caused primarily by thermal tails in the observations. Objects with 
thermal tails are labeled with a footnote and will be addressed in a 
future paper. 

We do not allow any past classifications to influence the classi- 
fications presented in this paper. We compare our classifications 
with the taxonomic classes reported by the European Asteroid Re- 
search Node (EARN) Near-Earth Asteroid Database 4 and find that 
most classifications from our survey are the same class or are within 
the same taxonomic complex as previously determined taxonomic 
classifications. Five objects have distinctly different classifications 
in our survey than appear in EARN. Three of these objects, (24761) 
Ahau, (154029) 2002 CY46, and (297418) 2000 SP43, have EARN 
classifications from broadband photometry observations (Hicks 
et al„ 2010; Somers et al., 2010; Ye, 2011). In these cases, the classi- 
fication differences are likely due to the low resolution of broadband 
photometry. We classify (297418) 2000 SP43 as V-type with our 
ExploreNEOs survey spectrum and with the MIT Joint Campaign 
spectrum, but the classification on EARN is Sr-type. This discrepancy 
is likely due to the fact that the previous classification comes from a 
visible wavelength spectrum with no near-infrared data. The final 
classification discrepancy is (138911) 2001 AE2. This object was pre- 
viously classified as T-type (Binzel et al., 2004b). The primitive, car- 
bonaceous-like, material implied by the T-type classification and the 
low delta-V orbit of the object made this object a candidate for 
spacecraft exploration. The classification for our observation is S- 
complex. The S-complex classification suggests the thermal history 
and evolution of the object is different from the primitive material 
suggested by the previous T-type classification. The T-type classifica- 
tion was determined with visible wavelength data and our near- 
infrared S-complex classification is unambiguous. This suggests that 
(138911) 2001 AE2 is not an appropriate target for missions to prim- 
itive objects. 

The spectrum for (433) Eros provides an example of a classifica- 
tion mismatch seen in our survey. Eros is often considered a quin- 
tessential S-type object; however the near-infrared spectra from 
the ExploreNEOs survey are formally classified in the Q-complex. 
Most of the MIT Joint Campaign observations are classified as S- 
complex. The key difference between the ExploreNEOs and MIT 
data is the availability of visible wavelength data. Our observations 
have no visible wavelength component included, while the MIT 
data uses visible wavelength data from Binzel et al. (2004c). We 
tested the robustness of these classifications by removing the vis- 
ible data from the MIT spectra and appending it to each of the 
ExploreNEOs observations. We reclassified each of these alternate 
spectra. All MIT spectra with no visible wavelength component 
classified as Q-complex (either Q/Sq or Q), while most ExploreN- 
EOs spectra with the visible wavelengths added classified as Sw 
(one classified as Sqw). This suggests that for near-infrared only 
observations, the Bus-DeMeo system returns more apparent Q- 
types than truly exist in the population. To address this problem 
quantitatively, visible spectra should be obtained for objects that 
have been classified as Q-complex using only near-infrared data. 
The classifications of the combined visible and near-infrared spec- 
tra can be compared against the classifications that use only near- 
infrared data. 


4 http://earn.dlr.de/nea/. 


2.4. Band parameters 

For objects belonging to the S-, Q-, and V-complexes, which 
show large 1- and 2-pm absorption features, certain spectral band 
parameters can be used to infer composition (e.g., Burns et al., 
1972; Adams, 1974; Gaffey et al., 1993; Dunn et al., 2010, 2013; 
Sanchez et al., 2012). These diagnostic absorption features are 
indicative of crystalline olivine and/or pyroxene. These minerals 
are common in various meteorites and asteroids. Olivine character- 
istically shows a broad absorption feature centered near 1.04- 
1.1 pm that is comprised of three overlapping bands, while pyrox- 
enes show two broad absorption features centered at 0.9-1 pm 
and 1.9-2 pm. These features are caused by Fe 2+ in the crystalline 
structure (e.g., Burns, 1993). For olivine-pyroxene mixtures, the 
wavelength of the 1 pm feature (Band I, BI) is a function of the rel- 
ative abundances and compositions of olivine and pyroxene, while 
the wavelength of the 2 pm feature (Band II, B1I) is a function of the 
pyroxene composition (e.g., Cloutis et al., 1986). The ratio of the 
area of Band II to Band I, or Band Area Ratio, is also a measure of 
the relative abundances of olivine and pyroxene. The depths of 
the bands can be used to investigate phase angle and space weath- 
ering effects under the assumption that mafic mineral abundance 
and particle size are the same (Sanchez et al., 2012; Reddy et al., 
2012c). The band parameters that work best for this analysis are 
1-pm Band Depth (BID), 1-pm Band Center (BIC), 2-pm Band Cen- 
ter (BIIC), and Band Area Ratio (BAR). For a very small number of 
taxonomically classified objects, the spectra are too noisy to calcu- 
late specific band parameters (e.g., (7358) Oze from 17 September 
201 1 does not have a BAR calculation). Calculated band parameters 
are shown in Tables 3 and 4. 

All spectral band parameters are calculated based on the proce- 
dures described in Gaffey et al. (2002). For consistency between 
our two reduction techniques, all spectra are normalized to unity 
at 1-pm. We find no change between the parameters calculated 
before and after the normalization. The spectra are displayed using 
this normalization in the appendix (Fig. 9). All parameters are cal- 
culated using a spectrum that is smoothed with a boxcar average of 
a chosen width. This ensures that no individual point can influence 
the calculation. Most objects are smoothed with a width of 10 
spectral points; however, some high signal-to-noise observations 
use a smaller value. To calculate band centers, a linear continuum 
between the two local maxima is calculated for each band. The lin- 
ear continua are divided out of each band and the center is calcu- 
lated by fitting a fourth-order polynomial to the bottom of the 
band and finding the local minimum. For some objects, the linear 
continuum that defines Band I will cut through the spectrum at 
slightly shorter wavelength than the local maximum at ~1.5 pm. 
This is defined as a “rollover” spectrum and in these cases a new 
local maximum is found from the continuum removed spectrum 
and the Band I Center calculation is repeated. For Band II, the red 
edge is nominally defined as the end of the spectrum. However, 
many spectra are noisy in this region and we remove some spectral 
points without changing the angle of the continuum in order to 
most accurately calculate the band centers and band areas. 

Band center errors are calculated via a Monte Carlo method. We 
generate offsets to each individual point in a spectrum by multi- 
plying the observational error by a normally distributed random 
number with a standard deviation of 1 . Each spectrum is altered 
20,000 times and a new band center is calculated for each trial. 
The band center error is defined as the standard deviation of these 
trials. We define a minimum band center error of 0.01 for all obser- 
vations to account for uncertainties that have not been quantified. 
Calculated uncertainties are only reported if they exceed this min- 
imum. Band II centers often have a much higher error than Band I 
centers. This is not only due to increased noise in Band II, but also 
to the shallowness and relative flatness of Band II, and therefore 
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Table 4 

Calculated band parameters for ExploreNEOs objects in the MIT-UH-IRTF Joint Campaign. As with Table 3, in order to distinguish between multiple observations of the same object the objects appear in the same order as they appear in 
Table 2. The table includes the same parameters as Table 3. 
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the polynomial fit, which allows for subtle polynomial shifts to re- 
sult in relatively large shifts in the band center calculations. 

Since some observations do not contain data that extend to 
~0.75 pm on the blue end, they do not contain a complete Band 
I. These Band I Centers are calculated using a linear continuum that 
begins at the start of the data, not at the local maximum blue-ward 
of Band I. Dividing Band I by this data-start continuum will pro- 
duce a differently sloped continuum-divided Band I than if it had 
been divided by an appropriate local maxima continuum. There- 
fore, there is an additional source of error on objects that do not in- 
clude data short ward of 0.75 pm. We quantify this error by 
removing the spectral data <0.82 pm for a few high signal-to-noise 
observations to simulate an observation taken with the SpeX di- 
chroic in the light path. We find that the BICs tend to shift to higher 
wavelength by values smaller than the defined minimum band 
center error of 0.01 . 

Band Area Ratio (BAR) is defined as the ratio of the area of Band 
II to the area of Band I. Each band area is calculated by using the 
trapezoidal rule to determine the area between the linear contin- 
uum and the continuum-divided spectral band. The error for each 
band area is calculated using the same Monte Carlo method used 
for band center errors. The error is defined as the standard devia- 
tion of the 20,000 altered area trials. The BAR error uses the indi- 
vidual area errors and propagates them with standard error 
analysis techniques. We define a minimum BAR error of 0.1 for 
all observations to account for uncertainties that have not been 
quantified such as the definition of the continuum slope and band 
depth changes due to phase angle. Calculated uncertainties are 
only reported if they exceed this minimum. This minimum error 
was also used by Reddy et al. (2010). The Band Area Ratio is not 
calculated for observations with no data <0.75 pm. Using a data- 
start continuum instead of a local maxima continuum removes a 
significant amount of area from the Band I area calculation. There- 
fore, the BARs calculated in this manner would not be directly 
comparable to the BARs calculated for asteroids using the complete 
Band I and the BARs calculated for laboratory meteorite spectra. 

Band depth is a measure of the ratio of the spectrum to the con- 
tinuum at the calculated band center. The depth of a band is diag- 
nostic of mafic mineral abundance and grain size, but is less 
diagnostic of composition than the other presented band parame- 
ters. The nominal band depth for an object can be modified by fac- 
tors such as phase angle and lunar-style space weathering. We 
present only Band I depth in this paper because the telluric absorp- 
tion remnants in Band II often make depth difficult to constrain. 
We use the smoothed, continuum-divided spectrum to determine 
Band I depth using the following equation: 


Band I Depth (%) 



* 100 


( 1 ) 


where R b is the reflectance spectrum value at the Band 1 center loca- 
tion and R c is the linear continuum value at the Band I center loca- 
tion (Clark and Roush, 1984). 


2.5. Temperature corrections of band parameters 

Surface temperature variations on asteroids shift the band cen- 
ters and change the width of the absorption feature (e.g., Moroz 
et al., 2000; Hinrichs and Lucey, 2002; Sanchez et al., 2012; Dunn 
et al., 2013). These temperature changes are of little importance 
for Main-Belt asteroids, but can be significant for the higher eccen- 
tricity near-Earth objects. Most importantly, in order to directly 
compare room temperature (~300 1<) laboratory spectra with 
observational spectra, temperature corrections are necessary. 

In order to apply temperature corrections to our sample, we 
need to identify the best laboratory matches to the two types of 


compositions in our band parameter sample. The S- and Q-complex 
objects are olivine-pyroxene assemblages. For these objects, we 
apply corrections derived from ordinary chondrites by Sanchez 
et al. (2012) and Dunn et al. (2013) using data from Moroz et al. 
(2000) and Hinrichs and Lucey (2002). The near-infrared spectra 
of the V-complex objects are dominated by the pyroxene features. 
For these objects, we apply corrections derived by Reddy et al. 
(2012c) using laboratory spectra of howardite, eucrite, and dioge- 
nite (HED) meteorites measured at various temperatures. All band 
parameters are corrected to room temperature values. Tables 3 and 
4 contain the temperature corrections and corrected band 
parameters. 

For each object, we calculate the average surface temperature 
(T) in Kelvin using the equation defined in Burbine et al. (2009). 

T = [(1 - A)L 0 /16>7£<T7ir 2 ] 1/4 (2) 

where A is the asteroid albedo, L 0 is the solar luminosity 
(3.827 x 10 26 W), r\ is the beaming parameter (in this case, assumed 
to be unity, Cohen et al., 1998), e is the asteroid’s infrared emissivity 
(assumed to be 0.9, Salisbury et al., 1991 ), a is the Stefan Boltzmann 
constant (5.67 x 10 s J s~ ] m~ 2 K~ 4 ), and r is the asteroid’s distance 
from the Sun in meters. When available, we use ExploreNEOs geo- 
metric albedos. Due to a small number of target list changes over 
the course of the Warm Spitzer observations, there are some objects 
on the ExploreNEOs target list for which albedos and diameters 
have not been measured. For those objects, we use the average al- 
bedo of the taxonomic complex (S = 0.26, Q_ = 0.29, and V = 0.42) de- 
rived by Thomas et al. (2011). All calculated corrections are added 
to the spectrally derived values. Corrections are on order of 0.001 
for BIC, 0.01 for BIIC, and 0.1 for BAR. We use Eq. (2) to calculate 
the asteroid surface temperatures to be consistent with previous 
asteroid spectroscopy temperature corrections (e.g. Burbine et al., 
2009; Moskovitz et al., 2010). The temperature calculated by this 
equation gives the temperature of a fast-rotating spherical asteroid 
with high thermal inertia and uniform temperature. This is known 
as the Fast-Rotating Model (FRM, Lebofsky and Spencer, 1989) 
and assumes that energy is absorbed over the side of the asteroid 
facing the Sun, but is radiated from the entire sphere. The temper- 
atures determined using Eq. (2) are lower than those calculated 
using the Standard Thermal Model (STM, Lebofsky et al., 1986). If 
we calculate the sub-solar temperature using the STM we get higher 
temperatures, which correspond to smaller corrections compared to 
the FRM temperature corrections. Both the FRM and STM correc- 
tions are on the order of our defined minimum band parameter er- 
rors and therefore the differences between them do not contribute 
significantly to the overall error budget. We note that the FRM is not 
consistent with asteroids in this size range with known rotation 
periods (e.g., Delbo et al., 2003; Muller et al., 2012; Emery et al., 
2010 ). 

The majority of our objects in the band parameter sample be- 
long to the olivine-pyroxene-rich S- and Q-complexes. Sanchez 
et al. (2012) shows that Band II is the most affected by temperature 
variations and that Band II center shifts to longer wavelength as 
temperature increases. The equation derived by Sanchez et al. 
(2012) for the temperature correction of the Band II center is the 
following: 

AB1IC (pm) = 0.06 - 0.0002 x T (K) (3) 

The calculated BIIC correction is added to the spectrally derived BIIC 
values. 

Sanchez et al. (2012) find the correction for Band I center to be 
very small in general (~0.003 pm) and do not calculate a temper- 
ature correction for it. Dunn et al. (2013) use different corrections 
from the LL chondrite Soko-Banja and the L chondrite Elenovka and 
also find the correction to be small. However, since the direction of 
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movement of the Band I center depends on the amount of olivine in 
the sample, the temperature corrections derived from these two 
ordinary chondrites do not shift in the same directions. Since both 
corrections are small, Dunn et al. (2013) apply a correction based 
on the average of the two corrections. We apply this method from 
Dunn et al. (2013) to our sample using the average of the following 
corrections: 

ABIC, LL correction (pm) = 0.000078 x T (K) - 0.02194 (4) 

ABIC, L correction (pm) = -0.000048 x T (K) + 0.01512 (5) 

We note that the constant in Eq. (4) is printed incorrectly as 0.21 94 in 
Dunn et al. (201 3). If we use 0.21 94 instead of 0.021 94, we calculate a 
temperature correction that is two orders of magnitude different 
from the corrections calculated using the complementary L chon- 
drites BIC correction (Eq. (5)). This would result in a correction of or- 
der 0.1 pm instead of <0.001 pm. Dunn et al. (2013) state that the 
corrections are consistent with Burbine et al. (2009) and are on order 
<0.001 pm. Therefore we use 0.02194 throughout our analysis. 

Sanchez et al. (2012) determine the correction for Band Area Ra- 
tio (BAR) to be the following: 

ABAR = 0.00075 x T (K) - 0.23 (6) 

Dunn et al. (2013) again use different corrections from the two 
meteorites and apply the average correction: 

ABAR, LL correction = 0.0009 x T (K) - 0.2666 (7) 

ABAR, L correction = 0.0005 x T (K) - 0.1708. (8) 

We find the Sanchez and Dunn BAR corrections to be nearly identi- 
cal and use the average of the two in our analysis. 

Band I depth was not addressed by Dunn et al. (2013) and San- 
chez et al. (2012) found no obvious correlation between Band I 
depth and temperature. No correction was applied. 

For the few pyroxene-rich V-types in our sample, we apply tem- 
perature corrections derived by Reddy et al. (2012c) using labora- 
tory spectra of howardite, eucrite, and diogenite (HED) meteorites 
measured at various temperatures. The band center corrections for 
howardites and eucrites were calculated using spectra from Hin- 
richs and Lucey (2002) and the corrections for diogenites were cal- 
culated using spectra from Schade and Wasch (1999). We apply an 
offset to each nominal V-type band center equal to the average of 
the two calculated corrections for the appropriate band. The cor- 
rections for Band 1 center are defined as the following: 

ABIC, Howardites k Eucrites (pm) 

= 0.01656-0.0000552 x T (K) (9) 

ABIC, Diogenites (pm) = 0.0000000017 x T 3 (K) 

- 0.0000012602 xT 2 (K) 

+ 0.0002664351 x T (K) 

-0.0124 (10) 

while the corrections for Band II center are the following: 

ABIIC, Howardites k Eucrites (pm) 

= 0.05067- 0.00017 xT(K) (11) 

ABIIC, Diogenites (pm) = 0.038544 - 0.000128 x T (K) (12) 

The resulting corrections are similar to those calculated with equa- 
tions derived by Burbine et al. (2009) who used spectra of enstatites 
at multiple temperatures measured by Moroz et al. (2000). Tempera- 
ture corrections for Band Area Ratio were not addressed in Reddy 
et al. (2012c), Burbine et al. (2009) or Moskovitz et al. (2010). There- 
fore, we do not correct the BARs of V-types for temperature variations. 


3. Results 

3.3. Band Parameters 

The results of our band parameter analysis are shown in Fig. 1. 
Panel (a) shows all band parameters for observed Q-, S-, and V- 
complex ExpIoreNEO targets in our survey (black circles) and in 
the MIT Joint Campaign for NEO Spectral Reconnaissance (red tri- 
angles). The Gaffey S-subtypes are shown in panel (b). These sub- 
types are the result of the S-complex containing a progression of 
mixtures that range from pure olivine (S(I)) to pure pyroxene 
(S(VII) and the basaltic acondrites (BA)). The comparison of the 
Gaffey S-subclasses and our calculated band parameters is shown 
in panel (c). Our data closely match the region defined by the con- 
tinuum of S-subtypes defined by Gaffey et al. (1993) and the oliv- 
ine-orthopyroxene mixing line determined by Clouds et al. (1986). 
However, it seems that a small offset to higher BAR could improve 
the fit. A calibration offset between laboratory derived data and 
telescopic data was proposed in Gaffey (1984) to explain a discrep- 
ancy between their results and calibrations in Adams (1974). It was 
not determined whether this correction was correct. We do not 
calculate an offset correction for our data, but do allow for an extra 
margin when investigating S(IV)-subclass members (panel (d), 
Section 3.5). 

Only one outlying object is excluded from this figure, (16834) 
1997 WU22 from sp77 of the MIT Joint Campaign. Visual inspec- 
tion of the spectrum shows an anomalous 2-micron band that is 
significantly different from other observations of 1997 WU22. 
We deem this spectrum to be anomalous and not indicative of 
the true composition and band parameters of the object. 


3.2. ExploreNEOs albedos and diameters 

We combine ExploreNEOs albedos with our calculated Band I 
Centers (BICs). For this analysis we use the albedos calculated via 
the method described in Mueller et al. (2011) and Trilling et al. 
(submitted for publication). The analysis produces an albedo and 
error bars generated by a Monte Carlo procedure. The calculations 
include a robust error analysis that considers the various sources of 
uncertainty present in the ExploreNEOs data and asymmetric error 
bars to capture the true uncertainty. The band parameter sample in 
this paper is taken from the Q-, S-, and V-complexes which are the 
higher albedo objects in our NEO sample. Due to the relatively 
large reflected light component for these objects, the albedo and 
diameter errors are correspondingly larger than for the low albedo 
objects (Harris et al., 2011; Mueller et al., 2011). 

To potentially identify correlations using the albedo data, we 
limit the sample to those with (T pV + and a pV - < 0.20. This error 
cutoff was chosen to eliminate the albedos with the highest error. 
Of the band parameters, we choose Band 1 Centers because they are 
the most accurate and precise measurements within our sample 
and the values can be used to determine the olivine and pyroxene 
chemistries (Dunn et al., 2010). 

As discussed in Trilling et al. (2010), the nominal albedo and 
diameter values for Eros are quite different from those reported 
by the NEAR mission and supporting ground-based observations. 
This discrepancy is due to the near pole-on geometry for the Spit- 
zer observations. Using a beaming parameter (»;) and absolute 
magnitude (H) that were more appropriate for the viewing geom- 
etry of the observation, Trilling et al. (2010) recalculated the albedo 
( p v ) and diameter (D) to be 0.23 and 22.3 km respectively. These 
recalculated values agree with previous observations of Eros (e.g. 
Harris and Davies, 1999; Li et al., 2004). We incorporate these 
recalculated values in our analysis. All other ExploreNEOs albedos 
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Fig. 1. (a) Band 1 Centers (BICs) versus Band Area Ratios (BARs) for all S-, Q-, and V-complex objects in our full spectral sample. Objects observed as part of the ExploreNEOs 
survey are displayed using black circles and those from the MIT Joint Campaign are shown using red triangles, (b) The Gaffey et al. (1993) S-subclasses and basaltic 
achondrites (BA) in B1C-BAR space. (c)The comparison of our spectral band parameters with the Gaffey et al. (1993) S-subclasses. The shape of the distribution of S-subclasses 
matches the distribution of objects in our analysis, (d) We use the S(IV) region as one method of identifying potential ordinary chondrites in our sample. This panel shows the 
S(1V) region with objects classified as S(IV) in our analysis. We include Eros and Itokawa because of their confirmed ordinary chondritic nature. (For interpretation of the 
references to color in this figure legend, the reader is referred to the web version of this article.) 


and diameters are the nominal values (Trilling et al., submitted for 
publication). 

One object, (141052) 2001 XR1, from the MIT Joint Campaign 
observations has multiple Warm Spitzer observations. This object 
was part of a subset of ExploreNEOs objects that were observed 
multiple times in order to study phase angle effects. We use the 
average of the four albedos and diameters in our analysis. We 
propagate the errors using standard error propagation techniques. 
The individual albedo and diameter values are shown in Table 4. 
The average values are p v = 0.22^ and D = 0.91 km. 

We add NEOWISE (Mainzer et al., 201 1 ) albedos and diameters 
to our analysis. For objects that lacked ExploreNEOs albedo and 
diameter solutions, we include the NEOWISE solutions. When the 
object appeared in both the ExploreNEOs and NEOWISE catalog, 
we average the albedos and diameters and propagated the errors 
using standard techniques in order to reduce the errors for these 
multiply observed objects. 

We find no correlation between our error-limited albedos or 
diameters (ExploreNEOs and NEOWISE) with respect to Band I Cen- 
ters (Fig. 2). Since object albedos can depend on object chemistry 
and olivine is more reflective on average than pyroxene (e.g. Reddy 
et al., 2011c, 2012a), one might expect to see a trend of increasing 
albedo with increasing Band 1 Center. We find no correlation be- 
tween albedo and BIC down to error limits of 0.05 in albedo. This 
suggests that processes such as space weathering and photometric 
effects due to various factors including regolith grain size have a 
non-negligible contribution to the observed albedo. 


1.05 - 


a 

3. 


0.95 


i 1 ~ — r 

This survey, E + N 
This survey, E only 
MIT Joint Campaign, E + N 
MIT Joint Campaign, E only 
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Fig. 2. We find no correlation between ExploreNEOs (E) and NEOWISE (N) albedos 
and Band I Center for our spectral sample. Objects observed as part of the 
ExploreNEOs survey are shown in black. The circles represent ExploreNEOs-only 
albedos and the squares are an average of ExploreNEOs and NEOWISE albedos. 
Objects from the MIT Joint Campaign are shown in red. The triangles represent 
ExploreNEOs-only albedos and the squares are an average of ExploreNEOs and 
NEOWISE albedos. The errors on the ExploreNEOs albedos of these objects are 
higher than the ExploreNEOs average since we are targeting the high albedo objects 
in this band parameter analysis. We place an upper limit on the errors o pV + and (j pV - 
of 0.20. We discard all observations with errors larger than this limit. We find no 
correlation down to error limits on o pV + and <r pV - of 0.05. (For interpretation of the 
references to color in this figure legend, the reader is referred to the web version of 
this article.) 
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3.3. Phase angle-Band Area Ratio correlations 

Band parameters can be influenced by factors not related to 
compositional variation. One important factor to consider is the 
phase angle of the observation. Phase angle (a) is the angle be- 
tween the light incident on the asteroid surface and the light 


(433) Eros 



(1036) Ganymed 



(1627) Ivar 



Fig. 3. Band Area Ratio (BAR) versus phase angle (a) for (433) Eros, (1036) 
Ganymed, and (1627) Ivar. For Eros and Ivar, BAR is anti-correlated with phase 
angle. Ganymed has an opposite correlation with a trend of increasing BAR with 
increasing phase angle. 


reflected from the surface towards an observer. Increased phase 
angle can result in an increase of spectral slope known as phase 
reddening and changes in the strength of the absorption bands. 

In Section 3.1 we note that a small offset to higher Band Area 
Ratio (BAR) could improve the match between our data and the 
laboratory-derived Gaffey S-subclasses and olivine-orthopyroxene 


(433) Eros 




(1627) Ivar 



Fig. 4. Near-infrared spectral slope versus phase angle (a) for (433) Eros, (1036) 
Ganymed, and (1627) Ivar. All three objects show evidence of phase reddening. The 
solid lines indicate the nominal correlations and the dashed lines indicate the 
correlations calculated without the outlying points. 
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mixing line. To investigate potential causes of this offset, we inves- 
tigate the effect of phase angle on BAR. We examine the 
relationship between phase angle and BAR on individual objects, 
not the full population, to understand the effects of varying compo- 
sitions and optical properties. Our data has not been corrected for 
phase angle effects although data presented in Sanchez et al. 
(2012) demonstrates that changes in phase angle will affect the 
spectral band depth. Ideally, we would correct BAR to a phase an- 
gle of 30° to best match laboratory spectra. Sanchez et al. (2012) 
did not find any correlations between increasing phase angle and 
BAR values. Our observations demonstrate that individual objects 
can show different spectral responses to changing phase angle 
and that there is not a universal phase angle-BAR correlation that 
can be applied to our spectral sample (Fig. 3). 

Three objects in our sample were observed at multiple phase 
angles: (433) Eros, (1036) Ganymed, and (1627) ivar. Eros was ob- 
served 1 1 times total by the ExploreNEOs and MIT Joint Campaign 
surveys. Using the phase angles (a) of the observations and the 
temperature corrected BARs, we calculate a least-squares linear 
fit of BAR = (-0.003 ± 0.002) x a + (0.445 ± 0.083) with R 2 = 0.180. 
R 2 is the coefficient of determination, which ranges from 0 to 1 
and indicates how well data fits a line (1 indicates the line perfectly 
fits the data). Ganymed also has 11 observations from the two sur- 
veys. We find a least-squares linear fit for Ganymed of 
BAR = (0.014 ± 0.005) x a + (0.773 ±0.131) with R 2 = 0.494. Ivar 
has 7 observations from the two surveys, which we use to find a 
least-squares linear fit of BAR = (-0.005 ± 0.004) x a + (0.454 ± 
0.098) with R 2 = 0.215. 

Eros and Ivar have similar phase angle-Band Area Ratio correla- 
tions: for these objects, BAR is anti correlated with phase angle. 
However, for Ganymed the correlation is opposite with a trend of 
increasing BAR with increasing phase angle. If most S-type objects 
behaved in a manner similar to Eros and Ivar, a correction from 
high phase angle to 30° could cause a shift in our BAR distribution 
to better match the Gaffey S-subclasses. Sanchez et al. (2012) saw 
similar correlations of BAR with respect to phase angle. Their two 
observations of Ganymed also displayed a correlation opposite to 
the other objects’. 

We find no consistent correlation between phase angle and BAR 
in our data. The large differences between the phase angle-BAR 
relationships for our three objects suggest that no individual cor- 
rection should be applied to all objects. Therefore, we do not cor- 
rect our sample for phase angle variations. 

3.4. Phase angle-near-infrared spectral slope correlations 

We also examine the relationship between phase angle and 
spectral slope in near-infrared wavelengths for (433) Eros, (1036) 
Ganymed, and (1627) Ivar. To investigate these correlations we 
use phase angles (a) and the near-infrared spectral slope between 
0.82 and 2.4 pm. We choose this wavelength range to be consistent 
for all spectra. Spectra from the ExploreNEOs survey do not include 
visible wavelength observations. Since we do not want to include 
data points from visible wavelength observations and several of 
the MIT Joint Campaign near-infrared spectra were taken with 
the SpeX dichroic in the light path, we select a lower wavelength 
of 0.82 pm. We choose the upper limit of wavelength to be 
2.4 pm to avoid contributions to the spectral slope from spectral 
noise at higher wavelengths. 

For the 1 1 observations of Eros we calculate a nominal least- 
squares linear fit of Near-Infrared Spectral Slope 
(NISS) = (0.005 ±0.001) x a + (0.180 ±0.053) with R 2 = 0.530. We 
find a least-squares linear fit for the 1 1 observations of Ganymed 
of NISS = (0.002 ± 0.001 ) x a + (0.096 ± 0.024) with R 2 = 0.302. We 
find a least-squares linear fit for the 7 observations of Ivar of 
NISS = (0.002 ± 0.002) x a + (0.231 ± 0.036) with R 2 = 0.186. 


Each of the three datasets has a single datapoint that lies far from 
the nominal correlation (Fig. 4). We also calculate each correlation 
without the single outlying datapoint. For Eros the outlying point 
is from the MIT Joint Campaign (spl05 on 28 December 2011). 
Without that observation we calculate a least-squares linear fit of 
NISS = (0.004 ± 0.001) x a + (0.197 ± 0.018) with R 2 = 0.867. For 
Ganymed the outlying point is from the ExploreNEOs survey on 17 
September 2011. Without the outlying observation we calculate a 
least-squares linear fit of NISS = (0.001 ± 0.001 ) x a + (0.1 07 ± 0.1 5) 
with R 2 = 0.099. The outlying point for Ivar is from the MIT Joint 
Campaign (sp77 on 27 January 2009). Without this outlying obser- 
vation we calculate a least-squares linear fit of 
NISS = (0.002 ± 0.001) x a(0.237± 0.019) with R 2 = 0.661. The R 2 
values improve for Eros and Ivar when the outlying points are 
removed. 

We find evidence for phase reddening for all three objects. Fig. 4 
shows the correlations for Eros, Ganymed, and Ivar. The solid lines 
indicate the nominal correlations and the dashed lines indicate the 
correlations calculated without the outlying points. 

3.5. Ordinary chondrite-Iike mineralogies 

A subset of the olivine-pyroxene mixtures represented in the 
Band I Center versus Band Area Ratio diagram represent mineralo- 
gies consistent with the ordinary chondrites. We identify the 
potential ordinary chondrites via two methods. First, we use the 
S(IV) region defined by Gaffey et al. (1993) to include the mineral- 
ogical assemblages typical of ordinary chondrites. Fig. Id shows 
the Gaffey S(IV) region and the S(IV) objects that were identified 
in our sample. This region, and similarly shaped regions, have been 
used by various other studies. For example, Dunn et al. (2010) used 
the S(IV) region defined by Gaffey et al. (1993) to determine the 
potential ordinary chondrites and de Leon et al. (2010) identified 
ordinary chondrite-like objects by calculating a similar, but slightly 
larger, region using the spectral band parameters of ordinary chon- 
drite spectra. We identify S(1V) objects by first selecting all objects 
with data and error bars wholly within the S(IV) subclass. Objects 
in the nearby S(III) region and with lower Band I Centers than the 
S(IV) region are included if the datum is close to the boundary and 
at least half of the associated error bar falls within the S(IV) region. 
We include some objects with higher Band I Centers than the S(IV) 
region if the datum and error fall wholly outside of all other Gaffey 
regions. 

We established these generous selection criteria to allow for the 
possible shift in BAR due to phase angle effects. All spectra of Eros 
and Itokawa that are not selected by this procedure are added to 
the sample due to the confirmed ordinary chondritic nature of 
these asteroids (e.g. Trombka et al., 2000; Nakamura et al„ 2011; 
Foley et al., 2006). Some of the Eros and Itokawa spectra fall in 
the S(II1) region to the left of the S(IV) region. For Eros this is con- 
sistent with the S(11I) and S(IV) spectra observed with the near- 
infrared spectrometer on the Near Earth Asteroid Rendezvous 
(NEAR) mission (Bell et al., 2002). 

The second method used to identify potential ordinary chon- 
drite analogues is by calculating the BIC-BAR regions that corre- 
spond to the Dunn et al. (2010) ordinary chondrite mineralogies. 
Dunn et al. (2010) used laboratory derived values of Fayalite (Fa, 
the iron-rich end member of the olivine solid solution series), Fer- 
rosilite (Fs, the iron-rich end member of the pyroxene solid solu- 
tion series), and oIivine/(olivine + pyroxene) in 48 equilibrated 
ordinary chondrites to define regions (with and without errors) 
for the H, L, and LL ordinary chondrites. The derived mineralogies 
and spectral band parameters for the same 48 ordinary chondrites 
enabled Dunn et al. (2010) to calculate the relationships between 
the mineralogical parameters (Fa, Fs, and ol/(oI + px)) and the spec- 
tral parameters (BIC and BAR). We use the defined H, L, and LL or- 
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Fig. 5. Our second method of identifying potential ordinary chondrites uses the 
Dunn ordinary chondrite regions calculated from the ranges of mineralogies 
determined in Dunn et al. (201 OJ. We calculated the regions in Band 1 Center-Band 
Area Ratio space and used the regions to determine potential ordinary chondrite 
types. Each ordinary chondrite region is labeled with its type. The Gaffey S(1I1) and 
S(IV) regions are included for comparison. The L and H regions closely approximate 
the S(1V) region, but the LL region is much larger than the corresponding subset of 
the S(IV) region. The Dunn LL region includes most of the Gaffey S(1II) regions. 

dinary chondrite regions to determine the corresponding region in 
B1C-BAR space. Objects are selected as potential ordinary chon- 
drites if the datum is inside one of the three regions or if at least 
half of the corresponding error bar enters one of the regions. 
Fig. 5 shows the three ordinary chondrite type regions compared 
to the S(IV) region and includes all the objects determined to be 
ordinary chondrite-like via this method. All of the objects inside 
the Gaffey S(1V) region are included as a subset of the Dunn 
region objects. Using this method, all observations of Eros and 
Itokawa are automatically included in the sample. The Dunn 
et al. (2010) H and L regions are reasonably consistent with the 
Gaffey S(IV) region. However, the LL region is much larger than 
the S(IV) region and includes the majority of the S(1II) region which 
Gaffey et al. (1993) suggests contains different meteorite analogs 
(urelites). 

We calculate the mineralogies for all of our potential ordinary 
chondrites identified via both the S(1V) region and the Dunn re- 
gions using the following three equations from Dunn et al. 
(2010): 

ol/(ol + px) = -0.242 x BAR + 0.728 (13) 

Fa (mol%) = -1284.9 x (BIC) 2 + 2656.5 x (BIC) - 1342.3 (14) 

Fs (moI%) = -879.1 x (BIC) 2 + 1824.9 x (BIC) - 921.7 (15) 

ol/(ol + px) is a modal ratio of olivine to olivine and pyroxene. Fa 

is the amount (in mol%) of fayalite and Fs is the amount (in 
mol%) of ferrosilite. The root mean square errors presented by 
Dunn et al. (2010) are 0.03 for ol/(ol + px), 1.3mol% for Fa, and 
1.4mol% for Fs. 

We classify each of our potential ordinary chondrite analogs 
into a subtype using the Dunn et al. (2010) ordinary chondrite min- 
eralogy regions with errors included (Fig. 5). The regions available 
for analysis are Fa (mol%) with respect to ol/(ol + px) and Fs (mol%) 
with respect to ol/(ol + px). We classify each object in both of these 
regimes to determine our final classification. Each object is deter- 
mined to be a type if it falls within a region that is unique to a spe- 
cific ordinary chondrite type. Objects that fall in the overlap 
regions are classified as potentially belonging to either class. In 
both regimes, there are several outliers near the H and LL ordinary 
chondrite regions. These are potentially not ordinary chondrite 
analogues and are classified as “H?" and “LL?”. These are the 



0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 


ol/(ol+px) 

Fig. 6. We use the Dunn et al. (2010) equations to derive the mol% of Fayalite (Fa) 
and Ferrosilite (Fs) and olivine/( olivine + pyroxene) ratio from band parameters of 
ordinary chondrite-like spectra. The nominal errors on the chemistry calculations 
are 0.03 for the ol/(ol + px) calculation, 1 .3 for the calculation of the mol% of Fa, and 
1.4 for the calculation of mol% of Fs. The top panel shows mol% Fayalite versus ol/ 
(ol + px). The bottom panel shows mol% of Ferrosilite versus ol/(ol + px). We 
determine a potential ordinary chondrite subtype (FI, L, LL) based on the location of 
the objects in each of these compositional regions. The quadratic equations impose 
an artificial upper limit on Fayalite and Ferrosilite within the LL ordinary chondrite 
subclass that can be clearly seen. 


objects whose BIC-BAR data placed them outside of the defined 
regions, but whose error bars did not allow for exclusion from 
the sample. In order to be conservative, the final classification gi- 
ven is the least certain of the two classifications. For example, if 
one calculation is H, but another is H/L, then the final classification 
is H/L. Fig. 6 displays the ordinary chondrite mineralogies of our se- 
lected objects. Objects in the Gaffey S(IV) sample are shown with 
black circles (this ExploreNEOs survey) or red ’ triangles (MIT Joint 
Campaign). As mentioned previously, all S(IV) objects are a subset of 
the Dunn regions sample. The Dunn classified objects that are not 
part of the S(IV) sample are displayed as squares and use the same 
color scheme to differentiate between this survey and the MIT Joint 
Campaign. The ordinary chondrite mineralogies calculated often 
identify Eros and Itokawa as LL ordinary chondrites. This suggests 
that, at least for the high signal-to-noise objects, this method is a 
reliable indicator of object composition. Objects with ordinary chon- 
drite-like spectra, their mineralogies, and ordinary chondrite sub- 
type analogs are listed in Tables 5 and 6. 


5 For interpretation of color in Fig. 6, the reader is referred to the web version of 
this article. 


236 


C.A. Thomas et al./ Icarus 228 (2014) 217-246 


Table 5 

Calculated mineralogies and ordinary chondrite types for all potential ordinary chondrite types in the ExploreNEOs spectral survey. Potential ordinary chondrites were 
determined using the S(IV) region and the Dunn et al. (2010) ordinary chondrite regions. Since the objects in S(IV) region are a subset of the Dunn classified objects, they are listed 
first. The table includes object number, object name, spectral type, corrected Band I Center (ABIC), corrected Band Area Ratio (ABAR), olivine/(olivine + pyroxene), Mol % of 
Fayalite (Fa), Mol % of Ferrosilite (Fs), ordinary chondrite type using Fa - ol/(ol + px), ordinary chondrite type using Fsol/(ol + px), and the final ordinary chondrite type 
classification using the least certain of the two classifications. 


Number 

Name 

Type 

ABIC 

ABAR 

ol/(ol + px) 

Fa 

Fs 

OC type Fa 

OC type Fs 

OC type 

S(1V) 

433 

Eros 

QJSq 

0.962 

0.345 

0.64 

24.2 

20.3 

L 

L 

L 

433 

Eros 

QJSq 

0.976 

0.507 

0.61 

26.5 

22.0 

L/LL 

L/LL 

L/LL 

433 

Eros 

Sq/Q 

0.983 

0.398 

0.63 

27.5 

22.7 

L/LL 

L/LL 

L/LL 

433 

Eros 

QJSq 

0.995 

0.324 

0.65 

28.8 

23.7 

LL 

LL 

LL 

1036 

Ganymed 

Sr/S 

0.919 

0.989 

0.49 

13.9 

13.0 

H? 

H 

H? 

1036 

Ganymed 

Sr 

0.929 

1.086 

0.47 

16.7 

15.0 

H 

H 

H 

1865 

Cerberus 

S 

0.930 

0.533 

0.60 

16.9 

15.1 

H 

H 

H 

1866 

Sisyphus 

S 

0.929 

0.549 

0.60 

16.6 

14.9 

H 

H 

H 

1866 

Sisyphus B 

S 

0.932 

0.991 

0.49 

17.6 

15.6 

H 

H 

H 

1866 

Sisyphus D 

S 

0.918 

1.131 

0.45 

13.7 

12.8 

H? 

H? 

H? 

1866 

Sisyphus F 

S 

0.918 

0.770 

0.54 

13.7 

12.8 

H? 

H? 

H? 

1917 

Cuyo 

Sr 

0.925 

0.733 

0.55 

15.5 

14.1 

H 

H 

H 

4283 

Cuno 

Sq/Q 

0.979 

0.554 

0.59 

26.9 

22.3 

L/LL 

L/LL 

L/LL 

5143 

Heracles 

a 

1.016 

0.297 

0.66 

30.4 

24.9 

LL 

LL 

LL 

5626 

1991 FE 

Sq 

0.943 

0.588 

0.59 

20.1 

17.4 

H 

H 

H 

5646 

1990 TR 

a 

1.025 

0.247 

0.67 

30.7 

25.2 

LL 

LL 

LL 

7358 

Oze 

S 

0.924 

1.114 

0.46 

15.3 

14.0 

H 

H 

H 

11066 

Sigurd 

Sr 

0.923 

1.105 

0.46 

15.0 

13.8 

H 

H 

H 

12711 

Tukmit 

a 

0.979 

0.307 

0.65 

26.9 

22.3 

LL 

LL 

LL 

12923 

Zephyr 

S 

0.942 

1.084 

0.47 

19.9 

17.3 

H 

H 

H 

16834 

1997 WU22 

Sq 

0.961 

0.476 

0.61 

23.9 

20.1 

L 

L 

L 

19764 

2000 NF5 

Sq 

0.971 

0.397 

0.63 

25.7 

21.4 

L/LL 

L/LL 

L/LL 

27346 

2000 DN8 

a 

1.028 

0.356 

0.64 

30.7 

25.3 

LL 

LL 

LL 

36284 

2000 DM8 

Sq 

1.018 

0.291 

0.66 

30.4 

25.0 

LL 

LL 

LL 

53789 

2000 EDI 04 

Sq 

0.949 

0.766 

0.54 

21.6 

18.5 

L 

H/L 

H/L 

68216 

2001 CV26 

Sq 

0.934 

0.536 

0.60 

18.1 

15.9 

H 

H 

H 

68216 

2001 CV26 

S/Sq 

0.949 

0.747 

0.55 

21.6 

18.4 

L 

H/L 

H/L 

68216 

2001 CV26 

Sq 

0.960 

0.427 

0.62 

23.8 

20.1 

L 

L 

L 

85839 

1998 Y04 

a 

1.001 

0.271 

0.66 

29.4 

24.2 

LL 

LL 

LL 

86067 

1999 RM28 

a 

0.966 

0.438 

0.62 

24.9 

20.9 

L/LL 

L/LL 

L/LL 

88254 

2001 FM129 

a 

0.976 

0.604 

0.58 

26.4 

22.0 

L/LL 

L/LL 

L/LL 

100926 

1998 MQ 

Q. 

0.945 

0.564 

0.59 

20.7 

17.8 

H 

H/L 

H/L 

137032 

1998 UOl 

Sq 

0.968 

0.485 

0.61 

25.1 

21.0 

L/LL 

L/LL 

L/LL 

137062 

1998 WM 

Sq 

0.930 

0.706 

0.56 

17.0 

15.2 

H 

H 

H 

138883 

2000 YL29 

a 

0.958 

0.385 

0.63 

23.3 

19.7 

L 

L 

L 

141498 

2002 EZ16 

Sq 

0.967 

0.512 

0.60 

25.0 

20.9 

L/LL 

L/LL 

L/LL 

143381 

2003 BC21 

Sr/S 

0.924 

0.937 

0.50 

15.3 

13.9 

H 

H 

H 

152931 

2000 EA107 

a 

1.010 

0.283 

0.66 

30.0 

24.7 

LL 

LL 

LL 

159402 

1999 AP10 

Sq 

0.950 

0.663 

0.57 

21.8 

18.6 

L 

H/L 

H/L 

163697 

2003 EF54 

Q 

0.965 

0.476 

0.61 

24.7 

20.7 

L/LL 

L/LL 

L/LL 

207945 

1991 JW 

a 

0.991 

0.323 

0.65 

28.4 

23.4 

LL 

LL 

LL 

214088 

2004 JN13 

Sq 

0.981 

0.566 

0.59 

27.1 

22.5 

L/LL 

L/LL 

L/LL 

218863 

2006 W0127 

Sq 

0.990 

0.278 

0.66 

28.3 

23.3 

LL 

LL 

LL 

219071 

1997 US9 

a 

1.024 

0.248 

0.67 

30.6 

25.2 

LL 

LL 

LL 

220124 

2002 TE66 B 

Sq 

0.957 

0.465 

0.62 

23.2 

19.6 

L 

L 

L 

265187 

2003 YS117 

Q 

0.972 

0.460 

0.62 

25.9 

21.6 

L/LL 

L/LL 

L/LL 

297418 

2000 SP43 

V 

0.936 

1.126 

0.46 

18.5 

16.2 

H 

H 

H 


2000 TJ1 

Sq 

0.946 

0.566 

0.59 

20.8 

17.9 

H 

H/L 

H/L 

Dunn regions 

1036 

Ganymed 

Sr 

0.923 

1.226 

0.43 

15.0 

13.8 

H 

H 

H 

1627 

Ivar 

S 

1.024 

0.616 

0.58 

30.6 

25.2 

LL 

LL 

LL 

1943 

Anteros 

Sq 

0.987 

0.631 

0.58 

28.0 

23.1 

LL 

L/LL 

L/LL 

3122 

Florence 

Sq/Q 

1.030 

0.499 

0.61 

30.7 

25.3 

LL 

LL 

LL 

6239 

Minos 

QJSq 

0.981 

0.267 

0.66 

27.2 

22.5 

LL 

LL 

LL 

8567 

1996 HW1 

a 

1.002 

0.536 

0.60 

29.5 

24.2 

LL 

LL 

LL 

23187 

2000 PN9 

Sq 

1.031 

0.646 

0.57 

30.8 

25.3 

LL 

LL 

LL 

39572 

1993 DQ1 

a 

0.985 

0.609 

0.58 

27.7 

22.9 

LL 

L/LL 

L/LL 

96590 

1998 XB 

Sq 

1.000 

0.564 

0.59 

29.3 

24.1 

LL 

LL 

LL 

137125 

1999 CT3 

a 

1.038 

0.680 

0.56 

30.7 

25.4 

LL 

LL 

LL 

141498 

2002 EZ16 

Sq 

0.990 

0.158 

0.69 

28.3 

23.3 

LL 

LL 

LL 

198856 

2005 LR3 

Sq 

0.995 

0.217 

0.68 

28.8 

23.7 

LL 

LL 

LL 

198856 

2005 LR3 

a 

1.005 

0.188 

0.68 

29.7 

24.4 

LL 

LL 

LL 

207945 

1991 JW 

a 

1.001 

0.143 

0.69 

29.4 

24.2 

LL 

LL 

LL 

220124 

2002 TE66 A 

Sq/Q 

0.976 

0.286 

0.66 

26.5 

22.0 

LL 

LL 

LL 

347813 

2002 NP1 

a 

0.961 

0.235 

0.67 

24.0 

20.2 

LL? 

LL? 

LL? 


2001 HC 

Sq 

0.968 

0.246 

0.67 

25.2 

21.1 

LL 

LL 

LL 
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Table 6 

Calculated mineralogies and ordinary chondrite types for all potential ordinary chondrite types for ExploreNEOs objects in the MITJoint Campaign. As with Table 5, the S(IV) 
objects are listed first. The table includes the same parameters as Table 5. 


Number 

Name 

Type 

ABIC 

ABAR 

ol/(ol + px) 

Fa 

Fs 

OC type Fa 

OC type Fs 

OC type 

S(1V) 

433 

Eros 

Sw 

0.950 

0.362 

0.64 

21.7 

18.5 

L 

L 

L 

433 

Eros 

Sw 

0.955 

0.378 

0.64 

22.7 

19.3 

L 

L 

L 

433 

Eros 

Sw 

0.960 

0.440 

0.62 

23.8 

20.0 

L 

L 

L 

433 

Eros 

Sw 

0.980 

0.276 

0.66 

27.1 

22.4 

LL 

LL 

LL 

433 

Eros 

Sw 

0.990 

0.179 

0.68 

28.3 

23.3 

LL 

LL 

LL 

433 

Eros 

Sw 

0.995 

0.339 

0.65 

28.8 

23.7 

LL 

LL 

LL 

433 

Eros 

Sqw 

1.020 

0.154 

0.69 

30.5 

25.1 

LL 

LL 

LL 

1036 

Ganymed 

S 

0.925 

0.690 

0.56 

15.6 

14.1 

H 

H 

H 

1036 

Ganymed 

S 

0.930 

0.869 

0.52 

16.9 

15.1 

H 

H 

H 

1627 

Ivar 

Sr 

0.949 

0.503 

0.61 

21.6 

18.5 

H/L 

H/L 

H/L 

1627 

Ivar 

S 

0.954 

0.497 

0.61 

22.6 

19.2 

L 

H/L 

H/L 

1916 

Boreas 

Sw 

0.955 

0.411 

0.63 

22.8 

19.3 

L 

L 

L 

1917 

Cuyo 

Sv 

0.940 

1.152 

0.45 

19.5 

16.9 

H 

H 

H 

3102 

Krok 

Sqw 

1.000 

0.465 

0.62 

29.3 

24.1 

LL 

LL 

LL 

5587 

1990 SB 

Sr 

0.925 

0.911 

0.51 

15.5 

14.1 

H 

H 

H 

5626 

1991 FE 

S 

0.930 

0.496 

0.61 

16.9 

15.1 

H 

H 

H 

6455 

1992 HE 

Srw 

0.935 

0.648 

0.57 

18.1 

16.0 

H 

H 

H 

6455 

1992 HE 

Sqw 

0.980 

0.300 

0.66 

27.0 

22.4 

LL 

LL 

LL 

11066 

Sigurd 

S 

0.925 

0.833 

0.53 

15.4 

14.1 

H 

H 

H 

11398 

1998 YP11 

Sr 

0.925 

0.603 

0.58 

15.7 

14.2 

H 

H 

H 

15745 

1991 PM5 

S 

0.945 

0.583 

0.59 

20.7 

17.8 

H 

H/L 

H/L 

16834 

1997 WU22 

S 

0.960 

0.406 

0.63 

23.7 

20.0 

L 

L 

L 

25143 

Itokawa 

Sqw 

0.970 

0.238 

0.67 

25.6 

21.3 

LL 

LL 

LL 

35107 

1991 VH 

Sq 

0.980 

0.311 

0.65 

27.1 

22.4 

LL 

LL 

LL 

36284 

2000 DM8 

Sq 

0.985 

0.396 

0.63 

27.7 

22.9 

LL 

L/LL 

L/LL 

53435 

1999 VM40 

Srw 

0.930 

1.085 

0.47 

17.0 

15.2 

H 

H 

H 

54690 

2001 EB 

S 

0.940 

0.509 

0.60 

19.5 

16.9 

H 

H 

H 

89355 

2001 VS78 

Sr 

0.925 

0.852 

0.52 

15.6 

14.2 

H 

H 

H 

137062 

1998 WM 

Sr 

0.915 

0.702 

0.56 

12.5 

12.0 

H? 

H? 

H? 

138258 

2000 GD2 

Sq 

1.010 

0.461 

0.62 

30.1 

24.7 

LL 

LL 

LL 

141052 

2001 XR1 

Sq 

0.960 

0.356 

0.64 

23.9 

20.1 

L 

L 

L 

143651 

2003 QO104 

Q 

0.996 

0.296 

0.66 

28.9 

23.8 

LL 

LL 

LL 

152931 

2000 EA107 

0. 

1.005 

0.274 

0.66 

29.7 

24.4 

LL 

LL 

LL 

200840 

2001 XN254 

S 

0.950 

0.695 

0.56 

21.9 

18.6 

L 

H/L 

H/L 

Dunn regions 

1627 

Ivar 

Sqw 

0.960 

0.213 

0.68 

23.8 

20.0 

LL? 

LL? 

LL? 

1943 

Anteros 

S 

1.010 

0.525 

0.60 

30.0 

24.7 

LL 

LL 

LL 

4183 

Cuno 

0. 

1.014 

0.179 

0.68 

30.3 

24.9 

LL 

LL 

LL 

5131 

1990 BG 

Sq 

1.018 

0.122 

0.70 

30.4 

25.0 

LL 

LL 

LL 

5143 

Heracles 

Q 

1.009 

0.148 

0.69 

30.0 

24.6 

LL 

LL 

LL 

6239 

Minos 

Sq 

0.985 

0.196 

0.68 

27.7 

22.9 

LL 

LL 

LL 

8567 

1996 HW1 

0. 

0.995 

0.202 

0.68 

28.8 

23.7 

LL 

LL 

LL 

35107 

1991 VH 

Sq 

1.005 

0.245 

0.67 

29.7 

24.4 

LL 

LL 

LL 

207945 

1991 JW 

Sq/Q 

0.970 

0.275 

0.66 

25.5 

21.3 

LL 

LL 

LL 

219071 

1997 US9 

0. 

1.030 

0.201 

0.68 

30.7 

25.3 

LL 

LL 

LL 


We have a total of 82 objects in the Gaffey S(IV) sample, which 
breaks down into 4.9% “H?”, 29.3% H, 11.0% H/L, 13.4% L, 14.6% L / 
LL, and 26.8% LL. The sample determined using the Dunn ordinary 
chondrite regions contains 109 objects and breaks down into 3.7% 
“H?”, 22.9% H, 8.3% H/L, 10.1% L, 12.8% L/LL, 40.4% LL, and 1.8% 
“LL?”. We discuss the sample statistics and the comparison to 
the meteorite population in Section 4. 

3.6. Spectral variation: Sisyphus and Cuyo 

As part of the survey, we obtained rotationally resolved spectra 
of two near-Earth Asteroids: (1866) Sisyphus and (1917) Cuyo. Sis- 
yphus was observed once on 1 May 2011, in six sequential obser- 
vations on 10 June 2011, and in three sequential observations on 
24 June 2011. Cuyo was observed in six sequential observations 
on 24 June 2011 and once on 14 August 2011. Fig. 7 shows Band 
I Center for Sisyphus (top panel) and Cuyo (bottom panel) for our 
sequences of observations. The x-axis is not directly correlated 
with time. Each set of observations are connected and presented 


sequentially. Additional observations for each object are off- 
set along the x-axis. 

The Band I Centers of the Sisyphus spectra vary from 0.915 to 
0.932 pm among the ten observations, implying that the amounts 
of fayalite and ferrosilite vary across the surface. While the 
variations are not greater than the errors associated with the 
BICs, it should be noted that the errors on BIC are artificially 
limited to be equal to or greater than 0.01. This error limit 
accounts for many unquantified sources of error, but this limit 
is much higher than the nominal Sisyphus BIC errors which are 
~0.001 -0.002. The nominal errors are low due to the high sig- 
nal-to-noise of the observations that also enable us to separate 
the observations into rotationally resolved components instead 
of requiring the co-addition of the frames over the rotating body. 
The low nominal errors suggest that Sisyphus likely contains 
compositional changes over its surface. The variation corresponds 
to the difference between the S(VII) region (for the 0.915 pm BIC) 
and the S(IV) region with H chondrite mineralogy (for the 
0.932 pm BIC). 
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(1866) Sisyphus 



(1917) Cuyo 



Fig. 7. Band I Center for Sisyphus (top panel) and Cuyo (bottom panel) for our 
sequences of observations. The x-axis is not directly correlated with time. Each set 
of observations is connected and presented sequentially. Additional observations 
for each object are offset in x. We see potential compositional variation on the 
surface of Sisyphus. Any variation is within the errors, but we set a hard minimum 
on the errors in order to better include several unquantified sources of error. The 
differences in BIC exceed the nominal calculated errors for the observations. Cuyo 
does not show indications of compositional variation on the surface. 

Cuyo shows a much smaller variation in Band 1 Center from 
0.922 to 0.926 pm. These BICs are also presented with the error 
limit of 0.01. Even within the low nominal Cuyo BIC errors 
(~0.001 -0.002, similar to those of Sisyphus), there is no evidence 
for compositional variation across the surface of Cuyo. 

4. Discussion 

We found no correlation between albedo and Band 1 Center. The 
absence of a correlation could be due to the large uncertainties in 
the ExploreNEOs albedo calculations (±50%, Harris et al., 2011). 
Other datasets should be used to confirm the lack of a relationship 
between these parameters before it is concluded that another pro- 
cess, such as space weathering or regolith evolution, is masking the 
nominal albedo of the grains. 

Our investigation of the relationship between phase angle and 
Band Area Ratio for (433) Eros, (1036) Ganymed, and (1627) Ivar 
confirms the phase angle-BAR correlations seen in Sanchez et al. 
(2012). The Sanchez et al. (2012) sample contained only two or 
three observations for each object and most correlations were flat 
or slightly negative in slope. Our results for Eros (11 spectra) and 
Ivar (7 spectra) show similarly negative slopes. Ganymed is an 
exception in both our ExploreNEOs sample and in the Sanchez 


et al. work with a positive slope that correlates increasing BAR with 
increasing phase angle. There is no a priori reason to expect the 
optical properties of Ganymed to be different than other objects. 
All three objects investigated in this work are classified as potential 
ordinary chondrites in our analysis. One difference of note that 
separates Ganymed from Eros and Ivar is the higher calculated 
BARs. This leads to Ganymed’s classification as a potential H chon- 
drite while Eros is often correctly identified as an LL chondrite (e.g. 
Trombka et al., 2000; McCoy et al., 2001) and Ivar has band param- 
eters associated with a range of potential ordinary chondrite types. 
While the BARs calculated for Ganymed are much larger than those 
of Eros and Ivar, Sanchez et al. (2012) presents band parameters of 
one other object at similar BAR, (4954) Eric, which shows a nearly 
flat, slightly negative correlation. This object also has band param- 
eters correlating to classification as a potential H chondrite. There- 
fore, it is unlikely that the differing slopes can be solely attributed 
to the compositions associated with the higher BARs calculated for 
Ganymed and Eric. Additional studies of the phase angle-BAR rela- 
tionship should focus on examining objects with a variety of BARs. 
The Sanchez et al. (2012) sample had 12 different objects and 10 of 
those had BARs less than or equal to ~0.6. Another possible expla- 
nation is that Ganymed is much larger than many other NEOs 
(D = 35.01 km, Usui et al., 2011). The variations observed in the 
BARs for Ganymed could be partially attributed to compositional 
heterogeneities on the surface. The potential for large objects to 
display surface heterogeneity does not apply to Eros, which dis- 
plays spectral homogeneity that was observed during the Near- 
Earth Asteroid Rendezvous (NEAR) mission (Bell et al., 2002). 

Our 11 near-infrared spectra of Ganymed suggest composi- 
tional heterogeneity on the surface. Five of the Ganymed spectra 
are classified as potential H chondrites. The other six fall either 
in the region with lower BIC than the S(IV) region (outside of any 
S subclass region), in the S(VI) region, or in the S(VII) region. Fie- 
ber-Beyer et al. (2011) classified Ganymed as S(VI) and suggested 
that it was a mesosiderite analog containing HED pyroxenes and 
NiFe metal. Our analysis does not confirm this result. If the entire 
range of mineralogies suggested by our band parameters exist on 
Ganymed, it is possible that the Fieber-Beyer et al. (201 1 ) spectrum 
observed only the S(VI)-type surface assemblages since their 
observations were taken in a time period (3.6 h) shorter than the 
rotational period of Ganymed (10.31 h, Skiff et al., 2012). 

The Dunn et al. (2010) equations for olivine and pyroxene 
chemistry used in Section 3.5 are quadratic. Unfortunately, the lo- 
cal maxima for both the fayalite and ferrosilite equations are with- 
in the range of laboratory derived chemistries of the LL ordinary 
chondrite. This implies that there is a range of fayalite and ferros- 
ilite values that exists in the ordinary chondrite population, but 
will never be predicted by the Dunn et al. (2010) equations. The 
maximum calculable value for fayalite is ~ 30.76 which corre- 
sponds to a Band I Center of 1.034 and the maximum calculable va- 
lue for ferrosilite is ~25.37 which corresponds to a Band I Center of 
1.038. This upper limit falls within the true LL ordinary chondrite 
region, not the larger region that includes the root square mean er- 
rors associated with spectrally derived values that is used in our 
analysis. The artificial limit caused by the turnover of the fayalite 
and ferrosilite equations is clear in Fig. 6 and does not appear to af- 
fect our ordinary chondrite class assignments. However, we note 
that the limit causes an artificial enhancement of objects with sim- 
ilar fayalite and ferrosilite values near the limit and an artificial 
dearth of objects with fayalite and ferrosilite values above the 
limit. 

One large unanswered question in the study of the composi- 
tions of near-Earth objects is the apparent mismatch between 
the relative proportions of the ordinary chondrite subclasses (H, 
L, LL) seen in meteorite falls compared to the relative proportions 
of the ordinary chondrite subclass analogues observed spectra- 
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scopically in near-Earth space. Burbine et al. (2002) calculate that 
the H, L, and LL ordinary chondrites are 34.1%, 38.0%, and 7.9% of 
all meteorite falls, respectively. This corresponds to H, L, and LL or- 
dinary chondrites being 42.6%, 47.5%, and 9.8% of ordinary chon- 
drite falls. However, the majority of near-Earth asteroids that are 
ordinary chondrite analogs have spectra consistent with LL chon- 
drite mineralogies (e.g. Vernazza et al„ 2008). The analysis of the 
47 NEOs with ordinary chondrite compositions by Dunn et al. 
(2013) found that H, L, and LL ordinary chondrite analogs are 
15%, 10%, and 60%, respectively, with objects that could not be dis- 
tinguished between L and LL accounting for 15%. de Leon et al. 
(2010) showed that the concentration of BAR and BIC values in 
the NEO population is consistent with a more significant presence 
of olivine in NEOs compared to Main Belt objects (see Fig. 5 of de 
Leon et al. (2010)). This concentration mostly occupies the LL or- 
dinary chondrite region. 

The 82 objects from our sample that fall in the Gaffey S(IV) sam- 
ple break down into 4.9% “H?”, 29.3% H, 11.0% H/L, 13.4% L, 14.6% 
L/LL, and 26.8% LL. If we include the “H?” objects in the H chondrite 
group, the Gaffey S(IV) sample contains 34.1% potential H chon- 
drites. The 109 objects from our sample in the Dunn ordinary 
chondrite regions sample break down into 3.7% “H?”, 22.9% H, 
8.3% H/L, 10.1% L, 12.8% L/LL, 40.4% LL, and 1.8% “LL?”. If we include 
the “H?” and “LL?” objects into their respective groups, the Dunn 
sample contains 26.6% potential H chondrites and 42.2% potential 
LL chondrites. The relative proportions of the ordinary chondrites 
in our samples (the S(IV) and Dunn regions) are between the H 
and L chondrite dominant meteorite statistics (Burbine et al., 
2002) and the LL chondrite dominant NEO statistics calculated by 
Dunn et al. (2013) and seen in Vernazza et al. (2008) and de 
Leon et al. (2010). Our samples have elevated percentages of po- 
tential LL ordinary chondrites compared to the meteorite fall sta- 
tistics, but do not have values as high as the Dunn et al. (2013) 
sample. Similarly, we do not show the same concentration of NEOs 
in BIC-BAR space (near the LL ordinary chondrites) as was seen in 
de Leon et al. (2010). Instead we show a broader distribution of 
NEOs in BIC-BAR space. This suggests that our objects are more 
distributed over the possible mineralogies than previous surveys. 

Our sample could be a step forward in finding the link between 
the meteorite population and previous NEO observations. How- 
ever, there are no clear physical distinctions between our sample 
and the Dunn et al. (2013) sample that would explain the discrep- 
ancies between them. There are 13 objects that are present in both 
our sample and the Dunn et al. (2013) sample. Most ordinary chon- 
drite classifications are consistent between the samples (there are 
two exceptions). One initial possibility was the sizes of the objects 
in each sample: small asteroids are the immediate precursors to 
meteorite falls and therefore should show mineralogies similar to 
meteorite falls. We can not directly compare diameters since not 
all of the objects in the Dunn et al. (2013) sample have thermal 
diameters. In order to be consistent, we compared the samples 
using H magnitudes from the Minor Planet Center as a proxy for 
diameter. The H magnitude ranges for both samples are nearly 
identical with a peak in the distribution around H ~ 16-17. Fig. 8 
shows the histograms with the ExpIoreNEOs sample in blue and 
the Dunn et al. (2013) sample in green. The ExpIoreNEOs sample 
is clearly not sampling a smaller size regime than the Dunn et al. 
(2013) sample. We also compared the orbital elements of the ob- 
jects in our survey to the objects in the Dunn et al. (2013) sample. 
There were no clear differences in semi-major axis (a), eccentricity 
(e), inclination (i), perihelion distance (q), or aphelion distance (Q). 

It is important to note that some observations of Eros and Itok- 
awa fall outside of the S(IV) region, but are included in the Dunn 
regions. Since these spectra have been observed with high signal- 
to-noise, this suggests that their positions outside of the S(IV) re- 
gion are not due to spectral errors and are an accurate representa- 



Fig. 8. We compared the distribution of H magnitudes (as a proxy for size) of our 
ExpIoreNEOs spectra to the spectral sample analyzed in Dunn et al. (2013) to 
determine if we were investigating different size regimes. We have different 
relative proportions of the ordinary chondrite types within our samples, which 
cannot be explained by differences in object sizes. 


tion of the Band 1 Centers, Band Area Ratios, and the resulting 
mineralogy. Additionally, our spectral samples have several objects 
with multiple spectral observations (e.g. Ganymed). In several 
cases, only a subset of the spectra for each object are classified as 
ordinary chondrite analogues using the S(IV) region. When the lar- 
ger Dunn region is considered, many of the spectra from these mul- 
tiply observed objects that did not fall into the S(IV) region are also 
classified as ordinary chondrite analogues. Ganymed, 2002 EZ16, 
and 2002 TE66 are among the objects whose full set of observa- 
tions were split between ordinary chondrite analogues and non-or- 
dinary chondrite-like objects prior to the inclusion of the Dunn 
regions. This suggests that classification using the larger Dunn re- 
gions (or at minimum, some subset of them that is larger than the 
S(IV) region) is a valid way to identify potential ordinary 
chondrites. 

It is clear that ordinary chondrite-like spectra can fall outside of 
the S(IV) region. However, it does not necessarily follow that the 
area outside the S(IV) region that is covered by the Dunn ordinary 
chondrite regions contains only ordinary chondrites. Gaffey et al. 
(1993) suggest that other mineralogies could be present in the re- 
gions labeled as S(III). The spectral signatures quantified by the 
band parameters used in this work are those caused by FeO. Other 
pyroxene and olivine chemistries that are not as well suited to 
band parameter analysis can overlap the ordinary chondrite re- 
gions in BIC-BAR space. For example, previous spectral studies 
(e.g. Gaffey et al., 1993; Sunshine et al., 2004) have cautioned that 
high-calcium pyroxene, a pyroxene variation that is not being mea- 
sured by the Dunn equations, can affect the calculated band area 
ratios. If this BIC-BAR overlap in compositions is primarily in the 
LL ordinary chondrite region then it could help explain the over 
abundance of NEOs with LL ordinaiy chondrite-like compositions 
compared to the meteorite fall population. 

5. Conclusions and future work 

We have completed a ground-based near-infrared spectral sur- 
vey of near-Earth objects in support of the ExpIoreNEOs Warm 
Spitzer project and complement our observations with near-infra- 
red spectra of ExpIoreNEOs targets from the MIT Joint Campaign 
for NEO Spectral Reconnaissance. We report taxonomic classifica- 
tions, spectral band parameters, and ordinary chondrite mineralo- 
gies for objects within our combined dataset of 340 observations of 
187 targets. We find no correlation between spectral band param- 
eters and ExpIoreNEOs albedos and diameters. This suggests that 
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various optical effects, such as space weathering, can change the 
nominal albedo expected from composition alone. Our analysis 
suggests that for spectra that contain near-infrared data but lack 
the visible wavelength region, the Bus-DeMeo system returns 
more apparent Q-types than exist in the NEO population. We find 
negative Band Area Ratio correlations with phase angle for (433) 
Eros and (1627) Ivar, but a positive Band Area Ratio correlation 
with phase angle for (1036) Ganymed, confirming the findings of 
Sanchez et al. (2012). We find evidence for spectral phase redden- 
ing for Eros, Ganymed, and Ivar. We identify ordinary chondrite 
subtypes for all potential ordinary chondrites in our sample using 
equations from Dunn et al. (2010). Our resulting proportions of H, 
L, and LL ordinary chondrites differ from those calculated for mete- 
orite falls and in previous studies of ordinary chondrite-like NEOs. 

One complementary investigation that has not been completed 
is a source region analysis for our various taxonomic types and or- 
dinary chondrite types. Future work will use the source region 
model described in Bottke et al. (2002) to do a full source region 
analysis of the sample. We will also incorporate our visible wave- 
length observations that were not described in this paper. 
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Appendix A 

We show our near-infrared spectroscopic observations from the 
SpeX component of our ExploreNEOs spectroscopy survey in Fig. 9. 
All spectra have been normalized to unity at 1-pm. 
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Fig. 9. This paper presents the near-infrared spectra observed at the NASA Infrared Telescope Facility (IRTF) using the SpeX instrument as part of the ExploreNEOs spectral 
survey. We include the spectra of all 125 observations of 92 ExploreNEOs targets. Each spectrum contains the object name and number, the UT date of observation, and, when 
appropriate, a letter designating the sequence order for all rotationally resolved observations. All spectra are normalized to unity at 1-pm. 
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